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INTRODUCTION 


In  past  years  aircraft  configuration  design  was  accomplished  mainly  by 
wind  tunnel  testing  while  flow-calculation  methods  contributed  little  because 
they  were  limited  to  simple  geometries  and  restricted  in  the  physical  processes 
that  they  represented.  This  virtually  exclusive  reliance  on  testing  had  dis¬ 
advantages  which  led  to  less-than-optimum  designs.  Each  potential  configura¬ 
tion  had  to  be  fabricated  as  a  wind  tunnel  model  with  corresponding  expense  and 
time  delay  and,  if  the  tests  suggested  design  changes,  the  process  should  have 
been  repeated,  though  on  occasions  this  could  not  be  done  because  of  the  time 
and  expense  involved  in  an  iteration  cycle.  Moreover,  tests  provide  incomplete 
information  in  that,  for  example,  static-measurements  may  be  restricted  to  a 
small  number  of  locations.  Force  and  moment  data  are  seldom  explained  in  terms 
of  flow  phenomena  and  the  extent  of  separated  regions  is  usually  not  deter¬ 
mined.  Finally,  there  is  the  necessary  scaling  from  model  to  full-scale 
vehicle  which  can  be  uncertain  so  that,  for  example,  separation  observed  on  a 
model  may  be  different  from  that  encountered  in  full  scale. 

Intensive  efforts  have  been  pursued  for  many  years  to  help  to  overcome  the 
limitations  of  calculation  methods  and  to  develop  them  to  represent  accurately 
the  flows  over  airplane  configurations.  Regardless  of  the  specific  algorithm 
employed,  the  cost  and  especially  the  elapsed  time  required  to  develop  a 
numerical  representation  of  a  given  configuration  is  much  less  than  that 
required  for  an  experimental  test,  once  such  a  basic  numerical  model  has  been 
developed,  many  variations  on  the  design  can  be  investigated  computationally 
and  only  the  most  promising  selected  for  wind-tunnel  testing.  A  flow  computa¬ 
tion  method  is  required  to  represent  the  geometry  of  the  airplane  and  the 
essential  properties  of  the  fluid  and,  to  achieve  the  objective,  limitations 
of  calculation  methods  have  been  removed  in  recent  years  by  the  development  of 
new  and  powerful  calculation  tools. 

There  are  three  possible  approaches  which  can  be  used  for  the  calculation 
of  the  viscous  flow  over  an  aircraft  configuration.  The  first  approach  makes 
use  of  the  Reynolds-averaged  Navier-stokes  equations  and  various  reduced  forms 
Including  the  so-called  parabolized  and  thin-layer  Navier-stokes  equations. 
Significant  advances  have  been  made  in  this  area,  for  example,  by  Shang  and 
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Scherr  [1]  who  made  the  first  attempt  to  numerically  simulate  the  flowfield 
around  a  complete  aircraft  by  solving  the  Navier-Stokes  equations.  To  demon¬ 
strate  the  feasibility  of  their  approach,  they  chose  the  hypersonic  research 

aircraft  X24C-10D  for  which  a  detailed  experimental  database  exists.  Using  a 

5 

mesh  system  around  5  x  10  nodes,  they  performed  impressive  calculations  at 
an  angle  of  attack  of  six  degrees  with  a  nominal  Mach  number  of  5.95,  and 
Indicated  areas  where  future  research  should  concentrate  to  make  this  approach 
more  efficient  and  practical. 

The  second  and  third  approaches  both  make  use  of  solutions  of  inviscid  and 
viscous  flow  equations  coupled  by  special  procedures.  The  second  approach  is 
based  on  the  two-dimensional  method  developed  by  Gilmer  and  Bristow  [2]  in 
which  an  empirical  inviscid  flow  model  is  used  to  represent  the  effects  of 
flow  separation  [2,3].  A  direct  boundary- layer  calculation  is  employed  up  to 
the  point  of  separation.  Downstream  of  this  point,  a  free  surface  is  intro¬ 
duced  to  model  the  separated  flow  region.  The  shape  and  the  length  of  the 
separation  zone  are  computed  by  satisfying  a  constant  pressure  boundary  condi¬ 
tion  on  the  surface  and  very  good  results  have  been  obtained  for  airfoils  at  a 
wide  range  of  angles  of  attack  including  stall. 

The  third  approach,  which  is  referred  to  as  the  interactive  boundary- layer 
approach,  uses  special  coupling  techniques  between  inviscid  and  viscous  flows 
and  novel  numerical  procedures.  The  particular  form  developed  by  Cebeci  et  al. 
[4]  is  very  general  and  allows  any  inviscid  flow  method  to  be  coupled  with 
solutions  of  the  boundary- layer  equations.  For  example,  in  its  application  to 
two-dimensional  subsonic  flows  over  airfoils,  it  employs  Halsey's  inviscid 
procedure  [5]  based  on  the  conformal  mapping  and  Fourier  analysis  techniques 
and  computes  the  flow  over  the  airfoil  and  wake.  Successive  viscous  sweeps  are 
performed,  after  each  of  which  the  external  inviscid  solution  is  recomputed, 
until  a  converged  solution  is  obtained.  The  boundary- layer  method,  which  is 
an  inverse  finite-difference  scheme  developed  by  Cebeci  [6],  uses  an  algebraic 
eddy-viscosity  formulation  due  to  Cebeci  and  Smith  [7]  and  is  able  to  compute 
flows  with  large  regions  of  separation  without  numerical  problems.  In  regions 
of  reverse  flow,  it  uses  the  FLARE  approximation  [8]  in  which  the  streamwise 
convective  terra  is  set  equal  to  zero  in  the  recirculating  region.  A  detailed 
description  of  the  method  and  of  its  application  to  a  range  of  airfoils  at 
angles  of  attack  up  to  and  including  stall  is  provided  in  [4],  The  results 


presented  in  Section  3.1  show  that  this  procedure  has  removed  a  major  obt> t a  ie 
in  that  flows  with  large  regions  of  separation  can  now  be  computed  accurately. 

A  comparison  between  the  interactive  procedure  [4]  and  a  method  based  on 
the  thin-Navier-Stokes  equation  has  been  reported  in  [9]  for  the  NACA  0012 
airfoil.  This  study  demonstrated  that  the  Navier-Stokes  approach  and  the 
interactive  boundary- layer  approach  gave  comparable  results  up  to  the  stall 
angle.  However,  the  Interactive  boundary- layer  approach  required  much  less 
computing  time. 

The  purpose  of  the  present  work  is  to  develop  a  general  method  for  comput 
ing  three-dimensional  flows  on  wings  with  leading-  and  trailing-edge  separa 
tion.  The  second  and  third  approaches,  described  above,  have  been  critically 
examined  to  determine  the  extent  to  which  they  can  fulfill  this  purpose.  It 
was  considered  that  the  second  approach  should  be  more  suited  to  large  regions 
of  trailing-edge  separation  whereas  the  third  approach  should  be  more  appro¬ 
priate  for  leading-edge  separation  bubbles.  Extensive  tests  showed  that  the 
second  approach  failed  to  converge  in  those  situations  for  which  it  was 
intended  and,  although  it  was  satisfactory  for  flows  with  small  trailing-edge 
separation,  considerable  effort  may  still  be  required  to  overcome  Its  limita¬ 
tions.  The  third  approach  on  the  other  hand,  has  proved  to  be  able  to  repre¬ 
sent  leading-  and  trailing-edge  separation  without  limitations,  as  discussed 
further  in  Section  3.0. 

The  remainder  of  this  report  has  been  arranged  in  five  sections.  In  F  - 
tion  2.0,  the  inviscid  method  is  briefly  described  prior  to  more  extensive 
descriptions  of  interactive  boundary- layer  methods  based  on  strip  theory  and 
quasi- three-dimensional  approximations.  The  viscous-flow  equations,  transfor¬ 
mations  and  solution  procedure  have  been  described  previously  in  [4]  for  two- 
dimensional  flows  and  in  [10]  for  quasi-three-dimensional  flows.  Results 
obtained  with  these  procedures  are  presented  in  Section  3.0  for  several  con¬ 
figurations  and  the  relative  merits  of  the  two  interactive  approaches  are  con¬ 
sidered.  The  computer  program,  which  embodies  both  methods  is  described  in 
Section  4.0  and  sample  input  data  for  the  three  contract  test  cases  are  pre¬ 
sented  in  Section  5.0.  The  report  ends  with  a  summary  of  the  more  important 
conclusions  in  Section  6.0. 


SECTION  II 

DESCRIPTION  OP  THE  INTERACTIVE  BOUNDARY  LAYER  METHODS 


The  method  described  here  combines  lrivlst Id  flow  and  Inverse  boundary 
layer  procedures  In  an  Interactive  manner  which  permits  the  calculation  of 
flows  with  leading  and  t r al 1 lng-edge  separation.  According  to  this  method, 
the  lnvlscld  method  of  Section  2.1  is  used  to  compute  the  flow  over  the  given 
configuration  with  a  zero  normal  velocity  boundary  condition.  The  resulting 
pressure  distribution  serves  as  a  boundary  condition  for  the  boundary  layer 
method  of  Section  2.2  which  computes  all  relevant  boundary  layer  parameters. 
Including  a  surface  blowing  distribution  to  simulate  the  displacement  thickness 
effect.  While  the  lnvlscldf low  calculations  (see  Section  2.1)  are  performed 
for  the  entire  configuration,  viscous  flow  calculations  are  presently 
restricted  to  the  wing  and  make  use  of  strip  theory  and  quasi  three  dimensional 
approxlmat Ions  discussed  In  Section  2.2. 

The  viscous  effects  computed  by  the  boundary- layer  method  are  then  used 
to  determine  a  distribution  of  normal  velocity  on  the  surface.  The  lnvlscld 
flow  method  Is  used  for  a  second  time  with  this  blowing  velocity  distribution 
as  a  boundary  condition  and  the  procedure  Is  repeated  until  convergence. 

2.1  lnvlscld  Plow  Method 

The  lnvlscld  flow  method,  which  is  the  first  order  surface  source  panel 
method  developed  by  Hess,  is  capable  of  computing  flow  about  completely 
arbitrary  configurations.  Because  of  Its  robustness  and  Its  availability,  It 
has  been  acquired  by  several  dozen  facilities  around  the  country  and  applied 
not  only  to  aircraft  but  also  to  ships,  submarines,  automobiles,  buildings,  and 
topographical  features.  The  three-dimensional  body  Is  represented  by  a  set  of 
plane  quadrilateral  panels,  as  shown  In  Fig.  1*.  A  three  dimensional  body  con 
slsts  of  lifting  and  nonllfting  portions.  A  lifting  portion,  such  as  a  wing 
or  pylon,  is  characterized  by  having  a  well-defined  trailing  edge,  from  which 
Issues  a  trailing  vortex  wake  and  along  which  the  Kutta  condition  Is  applied. 
Nonllfting  portions  of  the  body  lack  such  trailing  edges.  Every  panel  has  a 
constant  value  of  source  strength  a.  Panels  of  lifting  portions  have.  In 
addition,  a  quadratically  varying  doublet  strength  p.  Values  of  source 
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density  on  the  panels  are  Independent  parameters  available  for  satisfying  the 
boundary  condition.  The  chordwise  dipole  variation  over  the  panels  of  a  lift¬ 
ing  strip  (Fig.  1)  is  assumed  in  a  form  that  leads  to  favorable  numerics,  so 
that  the  dipole  strength  is  reduced  to  a  single  adjustable  parameter  -  the 
circulation  on  that  strip.  Thus  the  set  of  independent  dipole  parameters 
equals  the  number  of  locations  where  the  Kutta  condition  is  enforced,  i.e.,  at 
the  trailing-edge  segment  of  each  strip. 

The  nature  of  the  Kutta  condition  adopted  by  the  many  panel  methods  which 
are  currently  available  varies  greatly.  The  condition  adopted  by  Hess  [11]  is 
the  physically  meaningful  condition  of  equal  upper  and  lower  surface  pressures 
at  the  trailing  edge.  All  other  methods  make  use  of  other  derived  conditions, 
e.g.  a  prescribed  flow  direction  a  short  distance  downstream,  which  do  not 
guarantee  a  pressure  match  at  the  trailing  edge.  For  instance,  Margason  et  al. 
[12]  showed  that  pressure  mismatches  of  up  to  half  of  the  freestream  dynamic 
pressure  could  occur  from  such  alternate  forms  of  the  Kutta  condition.  Figure 
2,  which  is  taken  from  [12]  shows  this  clearly.  Since  we  are  interested  in  the 
computation  of  flows  for  which  the  behavior  of  the  boundary  layer  at  the  trail¬ 
ing  edge  can  have  a  significant  effect  on  the  overall  flow  solution,  it  is 
believed  that  the  approach  adopted  here  is  more  realistic. 


The  computational  procedure  is  as  follows.  On  each  panel  a  control  point 
is  selected  where  the  normal-velocity  boundary  condition  is  applied  and  where 
velocity  and  pressure  are  eventually  calculated.  Using  the  well-known  point- 
source  potential,  the  velocity  Induced  by  the  source  density  on  a  panel  at  a 
point  in  space  is 


v(panel)  =  [  JJ  grad(l/r)dS]ff (panel) 
panel 


(2.1) 


where  the  constant  value  of  a  has  been  taken  outside  the  integral  whose  value 

therefore  depends  only  on  geometry.  Suppose  N  panels  are  used  to  define  the 

body,  and  let  a ^  denote  the  value  of  source  density  on  the  jth  panel. 

The  velocity  induced  by  the  jth  panel  at  the  control  point  of  the  ith  panel  is 
♦  -* 

defined  as  v^jaj'  where  an  integral  of  the  form  of  Eq.  (2.1).  The 

integral  may  be  carried  out  analytically.  Explicit  formulas  are  contained  in 
[11],  but  they  are  too  lengthy  for  inclusion  here.  The  corresponding  normal 
velocity  at  the  ith  control  point  is 


*lj”)  •  "t  •  Jlj') 
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where  Is  the  unit  normal  vector  to  the  1th  panel.  Applying  the  normal 
velocity  boundary  condition  at  all  control  points  then  yields  the  following 
set  of  linear  algebraic  equations: 


N 

I 

J-l 


(  2  .  i  ) 


This  Is  the  numerical  approximation  of  the  integral  equation  that  expresses 
the  zero  normal-velocity  boundary  condition.  Once  the  have  been  deter 
mined,  the  disturbance  velocities  at  the  control  points  due  to  the  body  are 
given  by 


-♦  N  _> 

v  =  l  V  o  1  =  1.  2,  - N  (2.4) 

J  =  1  3  3 

While  the  addition  of  lift  In  two  dimensions  causes  no  significant 
Increase  in  the  complexity  of  the  problem,  the  problem  of  three  dimensional 
lifting  flow  Is  not  only  considerably  more  complicated  than  nonlifting  flow 
but  requires  assumptions  that  are  somewhat  arbitrary  [11].  The  main  features 
are  Illustrated  in  Fig.  1.  As  mentioned  above,  wings  or  other  lifting  portions 
of  the  configuration  are  characterized  by  having  trailing  edges  from  which 
Issue  trailing  vortex  wakes.  So-called  bound  circulation  Is  hypothesized  to 
lie  on  or  within  the  wing  surface,  with  strength  varying  In  both  the  chordwise 
direction  and  In  the  direction  parallel  to  the  trailing  edge  ("spanwlse"  In 
aerodynamic  Jargon).  The  variation  of  the  bound  circulation  In  the  chordwise 
direction  Is  predetermined,  while  the  spanwise  variation  is  adjusted  to  satisfy 
a  condition  of  smooth  flow  off  the  entire  trailing  edge,  and  the  trailing  vor 
ticlty  has  constant  strength  In  the  stream  direction  and  an  Initial  strength 
at  the  trailing  edge  equal  to  the  local  "spanwise"  derivative  of  bound  clrcula 
tlon  strength.  The  location  of  the  wake  once  it  leaves  the  trailing  edge  Is 
Initially  unknown,  which  Introduces  a  nonlinearity  Into  the  problem,  it  Is 
customary  simply  to  assume  a  wake  location,  because  In  most  problems  the  calcu 
lated  results  are  not  sensitive  to  the  details  of  the  wake  shape.  However,  the 
location  can  be  determined  by  iterating  the  calculation  If  necessary. 


With  regard  to  computing  effort  and  cost,  the  two  principal  tasks  assocl 

^ j  and  t  he 


2  ♦ 

ated  with  this  method  are  calculation  of  the  N  velocities  v 
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solution  of  the  get  of  linear  equations.  Bq.  (2.3).  Other  portions  of  the 
calculation  require  camper  at lwe ly  neqltqlble  comput Inq  efforts.  In  three 
dimensional  applications  r a t he r  larqe  numbeis  of  panels  are  employed.  Panel 


nuahers  over  3000  are  cc 


These  are  the  number  of  unknotats  actually 


solved  for  after  all  syaaetrles.  etc.,  have  been  utilized.  Accordingly .  It  Is 
econoalcally  Important  that  all  possible  efficiencies  be  eaployed  In  carrylnq 
out  the  two  principal  tasks.  Por  the  first  task  this  Is  accoaplished  by  uslnq 
staple  approximate  expressions  for  V  when  computing  the  Influence  of  distant 
panels  [13].  The  second  task  Is  speeded  up  by  the  use  of  block  iterative  mat 
rix  solutions,  as  developed  by  Clark  [14).  The  source  method  lends  Itself 
particularly  well  to  these  efficiencies,  without  which  no  three  dimensional 
method  Is  practical. 


In  updating  the  external  lnvlscld  velocity  to  account  for  the  boundary 
layer  displacement  effects,  two  approaches  can  be  considered.  The  first,  known 
as  the  displacement  surface  approach,  involves  modifying  the  geometry  to 
a  count  for  the  boundary  layer  displacement  thickness.  The  second  approach  Is 
to  introduce  a  blowing  velocity  on  the  surface  of  the  original  body  so  that 
the  dividing  stream  surface  of  the  lnvlscld  flow  approximates  the  displacement 
thickness  computed  by  the  boundary-  layer  procedure.  While  both  approaches 
nave  been  applied  In  our  studies,  the  first  approach  has  been  found  to  be  less 
suitable  because  It  Involves  the  computation  of  the  lnvlscld  flow  past  a  body 
with  a  finite  opening  at  the  trailing  edge.  Por  separated  flow,  the  size  of 
this  open  trailing  edge  can  become  significant  and  the  computation  of  the 
lnvlscld  flow  past  such  a  body  is  nonunique  unless  some  additional  conditions 
are  Imposed  along  this  trailing- edge  opening.  By  adopting  the  surface  blowing 
approach,  the  equivalent  trailing  edge  thickness  of  the  dividing  stream  surf¬ 
ace  Is  controlled  by  the  blowing  velocity  introduced  on  the  wing  surface  In  a 
physically  realistic  way  without  resorting  to  additional  boundary  conditions. 

The  blowing  velocity  computed  by  the  boundary- layer  procedure  is  applied 
on  the  wing  surface.  However,  as  pointed  out  above,  the  aim  is  to  model  the 
flow  past  the  displacement  surface.  Por  flows  Involving  large  traillng-edge 
separation,  it  has  been  shown  in  [4]  that  it  is  Important  to  both  evaluate  the 
external  lnvlscld  pressure  distribution  and  apply  the  Kutta  condition  on  the 
displacement  surface.  in  the  present  method  this  is  achieved  by  introducing 
additional  off  body  points  corresponding  to  each  control  point  for  which  the 


■j  vS  v  v_V  / 


boundary  layer  dtstrlbutton  has  been  computed.  These  off  body  points  are  used 
both  to  compute  the  final  pressure  distribution  and  to  apply  the  Kutta 
condit Ion. 


2.2  interact lvt  Viscous -Flow  Method 

An  accurate  prediction  of  the  flowfleld  over  a  wing  requires  the  calcula 
t Ion  of  flow  on  the  surface  and  In  the  wake.  The  airfoil  studies  conducted  In 
[4]  Indicate  that  the  wake  Influence  is  negligible  for  airfoils  at  low  angles 
of  attack  and  that  It  Is  sufficient  to  calculate  the  flow  on  the  airfoil  only, 
in  the  case  of  the  NACA  0012  airfoil,  for  example,  the  wake  effect  begins  to 
become  Important  for  angles  of  Incidence,  a,  greater  than  10* ,  and  Is 
Important  at  angles  of  incidence  approaching  stall,  which  for  this  airfoil  Is 
around  16°  as  shown  In  Fig.  3.  The  Inclusion  of  the  wake  effect,  which  Is  not 
considered  here,  reduces  the  flow  separation  on  the  surface  and  allows  the 
numerical  calculations  to  be  performed  with  less  difficulty.  In  the  present 
study,  we  Investigate  the  ability  of  the  Interactive  boundary- layer  scheme  to 
calculate  flows  with  massive  separation  and  perform  calculations  for 
configurations  at  high  angles  of  attack,  as  discussed  In  Section  3.0. 


2.2.1  Boundary-Layer  Squat  ions 

The  full  three-dimensional  boundary- layer  equations  and  their  boundary 
conditions  may  be  written  as  the  following  nonorthogonal  curvilinear 
coordinate  system  [10]: 


Continuity  Equation 


ddd 
.  (uh_  sin0)  i  T~  (wh.  sin8)  +  ~  (vh.h„  sln8)  =  0 
8x  2  3z  1  3y  1  2 


x-Moraentum  Equation 


u  3u 
h^  3x 


w  3u  3u  _ 

r~  t~  +  v  —  -  cot 9 

h  3z  3y 


2  a  2 
K^u  +  csc8  K2w  +  K  uw 


2 

esc  8  3p  cot8  csc8  3p  3  3u 
ph^  3x  ph2  37  3y  v  3y 


u’V  ) 


(2.5a) 


(2.5b) 
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z  Women turn  Equation 


u_  dw 
h^  3x 


W  dw 

+  —  —  ♦  V 

h2  dz 

dw  .  2 
~  +  csc8  K.u 
dy  1 

.  2 

-  cot8  K  w  +  K  , uw 

2  21 

cot 8  csc8 

Phl 

.  2«  , 

dD  CSC  8  dp 

dx  ph2  dz  + 

8  dw  — 

*  <v  VW) 

(2.5c) 

y  3  0: 

u,  v,  w  *  0 

(2.6a) 

y  *  4: 

u  *  u  (x,z). 

w  =  w  (x,z) 

(2.6b) 

Here  x,z  denote  the  coordinate  system  on  the  surface  of  the  body  and  y  is  the 
actual  distance  measured  normal  to  the  surface.  The  boundary- layer  equations 
and  boundary  conditions  for  this  system  according  to  first-order  boundary- layer 
theory  are  based  on  the  assumption  that  the  pressure  is  constant  across  the 
shear  layer  and  stress  gradients  in  directions  parallel  to  the  surface  are 
negligible  compared  with  those  normal  to  the  surface.  In  the  above  equations 
h^  and  h^  denote  the  metric  coefficients  and  6  denotes  the  angle  between 
the  coordinate  lines  and,  as  a  result  of  first-order  boundary- layer  theory, 
they  are  functions  of  the  surface  coordinates  x  and  z  only.  They  can  be 
obtained  from  the  relations  which  define  the  three-dimensional  body  in  the 
Cartesian  coordinate  system  x.y.z  by 


F(x,y,z)  =  0  (2.7) 

and  those  that  define  the  curvilinear  coordinate  system.  Thus  the  metric 
coefficients  and  the  angle  between  8  are  given  by 


-  a"  2  a"  2  a“  2 

hl  *  W  +  (dt)  +  (di> 

(2.8a) 

-  2  -  2  -  2 

A  -  <H>  *  «£»  *  <ff> 

(2.8b) 

0  (dx/dx) (dx/dz)  +  (dy/dx) (dy/dz)  +  (dz/dx)(dz/dz) 

h.h. 

(2.9) 

The  parameters  and  are  known  as  the  geodesic  curvatures  of  the 
curves  z  *  constant  and  x  =  constant,  respectively,  see  Pig.  4,  and  are  given 
by 
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(2.10a) 


i  a  3hi 

K1  =  hth2  sine  [ai  <h2  cos6>  -  3z  1 


1  d 

K2  *  h1h2  sln8  [3i  (hl  cos9)  “  aT1 


(2.10b) 


The  parameters  K^2  and  K21  are  defined  by 

K 


1  1  36  1  36 

12  '  ilSe  c'  Ki  -  ^  a!  *  cose  <k2  *  r  a?  1 


(2.11a) 


1  1  36  l  36 

C21  =  sin6  t_  K2  '  h2  3z  +  COS0  (K1  +  hL  3x)  1 


(2.11b) 


The  magnitude  of  the  velocity  vector  ut  in  the  boundary  layer  is  given  by 

2  2  1/2 

ufc  =  (u  +  w  +  2uw  cos6)  (2.12) 


In  the  present  study  we  use  two  reduced  forms  of  the  above  boundary- layer 
equations.  The  first  is  referred  to  as  the  quasi-three-dimensional 
boundary- layer  equations  in  which  the  flow  variations  with  respect  to  z  are 
neglected  so  that  the  equations  become 


(uh2  sin8)  +  (vhjh2  sin6)  *  0 


(2.13a) 


u  du  ^  3u 
hx  3x  v  3y 


..2  . «...  2  «...  esc  6  dp  .  3  ,  3u 

Kl“  cote  *  k2»  csc6  *  K12uw  .  -  —  ^  ^  <»  ^ 


U '  V  '  ) 

(2.13b) 


^  I?  ♦  V  I?  *  K.u2  csce  -  K„w2  cote  *  K..UW  -  cs=9  cote  I?  ♦  I-  (V  I?  -  ^ 


hx  3x  w  3y  “'1 


21v 


V 


3x  3y  3y 


) 


(2.13c) 


The  second  is  referred  to  as  the  strip-theory  approximation  which  essentially 
solves  the  well-known  two-dimensional  boundary- layer  equations. 


3u  3v 
3x  +  3y 


0 


u 


3u 

3x 


+  v 


3u 

ay 


du 
_ € 

ae  dx 


3_ 

ay 


(V 


du 

ay 


-  U  ’  V  '  ) 


(2.14) 

(2.15) 


In  this  approach,  the  strearawise  external  velocity  u  is  replaced  by  the 

e 

total  velocity  V  defined  by 
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(2.  It) 


.  2  2  _  1/2 
V  *  (u  *  w  +  2u  w  cos6) 
e  e  e  e 

and  Is  assigned  to  each  spanwise  strip. 

2.2.2  Interactive  Schama 

It  Is  well  known  that  the  boundary- layer  equations  are  singular  at 
separation  when  solved  for  a  prescribed  external  vel<  city  distribution.  They 
are  not  singular  at  separation,  however,  when  the  external  velocity  Is  computed 
as  part  of  the  solution  by,  for  example,  prescribing  a  displacement  thickness: 
this  Is  known  as  the  Inverse  problem  and  leads  to  solution  of  the  boundary- 
layer  equations  with  separation. 


For  flows  with  separation,  the  presence  of  the  upstream  velocity  component 
introduces  a  numerical  instability  into  the  solution  of  the  boundary- layer 
equations.  Over  the  years,  several  approaches  have  been  proposed  to  overcome 
this  problem.  One  popular  approach,  known  as  the  FLARE  approximation  after 
the  originators,  Flugge-Lotz  and  Revnher  [8],  neglects  the  longitudinal  con¬ 
vection  terms  u(du/dx)  and  u(dw/dx)  in  the  momentum  equations  (2.5). 

This  approximation  is  satisfactory  provided  that  the  size  of  the  separated 
region  remains  small.  However,  as  the  size  of  the  recirculation  region 
Increases,  this  approximation  becomes  less  accurate  and  requires  the  develop 
ment  of  additional  procedures  to  account  for  the  neglected  convective  terms. 

One  successful  scheme,  which  has  been  applied  to  two-dimensional  flows,  is 


referred  to  as  the  DUIT  procedure  [15,  16]  (Downstream,  Upstream  Iteration) 
and  requires  several  sweeps  through  the  recirculation  region,  with  this 
scheme,  the  FLARE  approach  is  used  to  compute  a  solution  within  the  reclrcula 
tlon  region,  and  then  in  successive  sweeps  the  u(du/dx)  term  is  progres 


sively  Introduced  until  it  is  fully  represented.  Another  successful  approach 
developed  by  Cebeci  [17]  makes  use  of  the  unsteady  boundary  layer  equations  In 
which  the  convective  term  is  introduced  into  the  equations  progressively  as  a 
function  of  time.  The  present  method  uses  only  the  FLARE  approach  for  both 
2D  and  quasi- 3-D  flows,  and  does  not  consider  the  effect  of  the  terms  neg 
lected  due  to  this  approximation. 

In  an  Interactive  boundary  Layer  scheme,  a  link  between  a  displacement 
thickness  and  external  fLow  is  provided,  and  two  types  of  procedures  have  been 
developed  for  this  purpose  for  two-dimensional  flows.  In  the  first  approach 
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[18-22],  the  solutions  of  the  boundary- layer  equations  are  computed  Initially 
for  a  prescribed  external  velocity  to  obtain  an  estimate  of  the  displacement 
thickness  A*(x)  distribution,  and  then  In  an  Inverse  mode  for  a  specified 
displacement  thickness  distribution  4*(x).  If  this  Initial  calculation 
encounters  separation,  A*(x)  Is  extrapolated  to  the  trailing  edge.  The  sub¬ 
sequent  boundary- layer  calculations  are  then  performed  In  an  Inverse  mode  to 
compute  the  blowing  velocity  needed  In  the  lnvlscld  flow  method,  in  general, 
this  procedure  leads  to  two  external  velocity  distributions,  u^(x)  derived 
from  the  Inverse  boundary- layer  solution  and  u^(x)  derived  from  the  updated 
approximation  to  the  lnvlscld  velocity  past  the  body  with  viscous  effects.  A 
relaxation  formula  In  the  form 


'V 


u  (x) 

A*V  1  ( x )  -  A*\x)  [1  ♦  w  (~~r~r 

el' 


1)]. 


0,  l.  2,  . 


(2.17) 


vv 

i 

lT- 


wheie  u>  denotes  a  relaxation  parameter.  Is  then  Introduced  to  define  an  updated 
displacement  thickness  distribution  and  to  obtain  new  solutions  of  the 
boundary- layer  equations  and  the  lnvlscld  flow  equations  so  that  the  Inter 
active  procedure  between  lnvlscld  and  viscous  flow  solutions  can  be  Iterated 
until  convergence  Is  achieved. 


rjw 

% 

iv:*’ 


Cvl- 
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The  second  approach  [23],  which  is  recommended  on  the  grounds  of  general¬ 
ity  and  physical  basis,  treats  the  external  velocity  u  (x)  and  the  dlsplace- 

e 

ment  thickness  A*(x)  as  unknown  quantities.  The  boundary- layer  equations 

are  solved  simultaneously  in  an  Inverse  mode  and  with  successive  sweeps  over 

the  body  surface.  For  each  sweep,  the  external  boundary  condition  Is  written 

o 

as  the  sum  of  the  lnvlscld  velocity  u  (x)  over  the  body,  and  a  perturba¬ 
tion  velocity  4u  (x),  that  Is, 


y  =  4,  u  (x)  =■  u  (x)  +  Au  (x) 

e  e  e 


(2.18) 


The  perturbation  velocity  Au  Is  computed  from  a  local  approximation 

based  on  a  thin  airfoil  theory  In  terms  of  the  local  blowing  velocity,  d/da 

(u  4*),  required  to  simulate  the  boundary  layer  thickness.  The 
e 

perturbation  velocity  is  written  as 


Au  (x)  =  ^  J  b  jjr  (u  4*)  - 

e  x  do  e  x  o 

x 

a 


(2.19) 
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where  the  Interaction  region  is  confined  to  [xft.  x^].  Introducing  a 
discrete  approximation  for  this  integral  enables  the  perturbation  velocity  to 
be  expressed  in  terms  of  the  geometric  coefficients  of  the  airfoil,  as 
discussed  in  [4]. 


This  two-dimensional  interactive  procedure  has  recently  been  extended 
under  an  AFOSR  contract  [10]  to  the  quasi- three-dimensional  equations  referred 
previously.  The  relationship  between  displacement  thickness  and  external 
velocity  needed  in  the  interactive  calculations  was  obtained  by  generalizing 
the  formulation  used  for  two-dimensional  flows.  The  irrotationality 
condition,  which  for  an  orthogonal  system  is 

h  [h2(we  +  4we)]  =  [fe  hl(ue  +  6ue)]  (2*20) 


was  used  to  provide  a  relationship  between  the  two  velocity  components  u 

e 

and  w^  and  shows  that  the  choice  of  computing  the  perturbation  velocities 

due  to  viscous  effects  is  not  arbitrary.  The  assumption  that  6u  (x)  is  a 

e 

function  of  x  alone  requires  that 


and  that 


T-  (6w  ) 

dx  e 


w 

e 


o 

w 

e 


(2.21) 


In  this  way  the  edge  boundary  conditions  for  a  quasl-three-dlmenslonal 
boundary- layer  flow  with  Interaction  are  given  by  Eqs.  (2.6)  and  (2.18). 


2.2.3  Turbulence  Model 

The  presence  of  Reynolds  shear  stress  terms  in  the  boundary- layer  equa¬ 
tions  requires  a  turbulence  model.  The  algebraic  eddy-viscosity  c  forrau- 

m 

latlon  of  Cebeci  and  Smith  [7]  is  used  here.  According  to  this  formulation 

for  two-dimensional  wall  boundary- layer  flows,  c  is  defined  by  two 

m 

separate  formulas,  given  by 


2  Idul 


{ 0 . 4y [ 1  -  exp(-y/A) ] }  Ytr  0  <  y  <  yc 


£m 


/  ( u  u)dy 

.  e 


I  Y 
'  tr ' 


y  c  .<  y  <  * 


(2.22a) 


(2.22b) 
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where 


A  =  26vut1, 
3u 

Tl  =  «*  37  ' 


.  1.1/2 
u  *  ( — )  , 

t  p  max 


1  +  5.5(y/S) 


(2.23) 


The  condition  used  to  define  y  is  the  continuity  of  the  eddy 

c 

viscosity;  from  the  wall  outward  Eq.  (2.22a)  is  applied  until  its  value  is 
equal  to  the  one  given  by  Eq.  (2.22b). 


In  Eq.  (2.22),  y  is  an  intermittency  factor  which  accounts  for  the 
transitional  region  that  exists  between  a  laminar  and  turbulent  flow.  It  is 
given  by 


Y 


tr 


1  -  exp[-G(x  -  xtr)  J  —  ] 

xtr  e 


(2.24a) 


Here  xtr  is  the  location  of  the  start  of  transition.  The  empirical  factor  G 
is 


G  » 


_1 _  JL  r'1-34 

1200  2  x„ 

v  tr 


(2.24b) 


where  the  transition  Reynolds  number  R  =  (u  x/v)  and  the  transition 

x  e  tr 

location  xtr  is  either  specified  or  cal£filated  from  the  empirical  formula  of 
Michel  [29] 


1.17«°-46  (1  *  ^S. 


•tr 


tr 


(2.24c) 


According  to  the  Cebeci  Smith  model,  the  parameter  a  in  Eq.  (2.22b)  is 
equal  to  0.0168  for  values  of  R-  greater  than  5000,  and  is  given  by  the 
expression  in  [24]  for  Rc  less  than  5000.  Studies  conducted  by  Head  and 

V 

his  associates  [25,26]  the  recent  experimental  data  of  Nakayama  [27]  and 
Simpson  et  al.  [28]  and  the  numerical  studies  of  Carter  [20]  indicate  that  in 
flows  with  strong  pressure  gradient,  the  value  of  a  should  also  be  changed 
when  Rg  5000.  Head  and  his  associates  recommend  that  a  in  Eq.  (2.22b) 
be  given  by 


a  =  a  F(r) 
eq 


(2.25a) 
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where 


a  =  0.002094  +  0.02672[1  -  exp(-0. 1163G) ] 
G  =  4.8285  (n  +  1.0717)1/2  +  1.8438 


F  = 


n 


6*  dp 
Tw  dx 


(5  -  4r>/(3  -  2r) 


r  <  1 


r  >  1 


(2.25b) 

(2.25c) 

(2.25d) 

(2.25e) 
(2.25f ) 


In  Eqs .  (2.25e)  and  (2.25f),  r  represents  the  ratio  of  the  local  rate  of  growth 
of  the  boundary  layer  to  the  rate  of  growth  of  the  corresponding  equilibrium 
layer. 


They  also  suggested  that  y  in  Eq.  (2.23)  be  replaced  by 

_ 2J) _ 

y  ~  1  -  erf [l/2(y/6  -  3)] 

where  8  Is  a  function  of  shape  factor  H. 


(2.25g) 


Simpson  et  al.  [28]  suggest  that 


o  =  0.0168/F2-5  (2.26a) 

Here  F  denotes  the  ratio  of  the  product  of  the  turbulent  energy  by  normal 
stresses  to  that  by  shear  stress  evaluated  at  the  location  where  shear  stress 
Is  maximum,  that  is. 


F  =  {(u'2  -  V2)  3u/3x  } 

-u'v'  (3u/3y)  (-u'v') 

max 


(2.26b) 


Before  Eq.  (2.26a)  can  be  used  in  Eq.  (2.26b),  an  additional  relationship  between 
2  2  -  - 

(u'  -  v*  )  and  (-u'v')  at  (-u'v')max  is  needed.  Here  we  assume  that  the 

ratio  in  Eq.  (2.26b), 


8 


-u'v' 


} 


(-u'v’) 


max 


(2.26c) 


is  a  function  of  R_  =  t  /(-u'v')  which,  according  to  the  data  of  Nakayama 
T  w  max  1 

[27],  can  be  represented  by 
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(2.26d) 


‘  1  ♦  21^(2  -  R.J.) 

for  Rt  <  1.0.  For  RT  >  1.0,  we  take  3  to  be 


(2.26e) 


Introducing  the  above  relationships  into  the  definition  of  F,  we  have  the 
following  expression  for  a,  according  to  Eg.  (2.26a), 


_ 0.0168 _ 

[1  -  3(3u/3x)/(3u/3y)]‘ 


(2.27) 


where  3  is  given  by  Eqs.  (2.26d)  and  (2.26e).  This  expression  is  used  here 
although  further  studies  are  clearly  required  to  evaluate  its  range  of 
validity.  Work  of  this  nature  is  in  progress. 

For  three-dimensional  flows,  the  above  formulation  was  generalized  as 

discussed  in  [16]  and  in  the  inner  region  e  is  defined  by 

m 


2r/3u 


duv  ,3w. 


(em>i  "  L  +  <fj>  +  2<p  <|f>  cos0] 


(2.28) 


where 


L  =  0.4y[l  -  exp(-y/A) ] ,  A  =  26  — 

ut 


T .  1/2 

/  tw. 

u,  •  <— > 


(2.29) 


Ttw  =  P  [(l7)  +  0  +  2(f^)  (a7>  cos0i 

1  w  J  w  1  w 


In  the  outer  region  em  is  defined  by 


(cra)  =  a  J  (ute  -  ut)dy 

O  0 


(2.30) 


and  a  is  given  by  Eq.  (2.27). 


2.2.4  Transformed  Equations 

The  equations  of  Section  2.2.1  may  be  solved  in  the  forms  presented  or 
expressed  in  other  forms  which  are  more  convenient  for  accurate  solution.  For 
two-dimensional  flows,  as  in  [4],  we  use  the  Falkner-Skan  transformation  in 
the  early  stages  of  the  flow 
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n  =  /u ft/vx  y,  ♦  3  /uevx  f(x,n) 


(2.31) 


With  primes  denoting  differentiation  with  respect  to  n  and  b  3  1  +  e^/v,  Eqs. 
(2.14)  and  (2.15)  and  their  boundary  conditions  can  be  written  In  the  form: 


(bf)’  +  ff-  +  m[l  -  (f*)2]  =  x(f*  -  f"  |f) 


n  =  0. 


n  =  n# 


f  =  f  3  0 


f  =  l 


(2.32) 

(2.33a) 

(2.33b) 


Here  f  Is  the  usual  definition  of  the  stream  function  that  satisfies  the 
continuity  equation. 


and  m  is  a  dimensionless  pressure-gradient  parameter, 

du 

x  _ e 

m  =  —  - — 
u  dx 
e 


(2.34) 


(2.35) 


This  transformation  provides  the  generation  of  initial  conditions  at  the 
stagnation  point  of  the  airfoil  and  allows  the  calculations  to  be  performed 
economically  and  accurately  around  the  leading  edge,  where  the  governing  equa¬ 
tions  are  being  solved  for  the  prescribed  external  velocity  distribution.  For 

Interactive  boundary- layer  calculations,  where  u  (x)  is  not  known,  a  constant 

e 

reference  velocity  u  Is  used  in  the  transformation 

o 


Y  3  /u  /vx  y,  ♦  =  /u  vx  F(x,Y) 
o  o 

In  terms  of  these  new  variables,  Eqs.  (2.14)  and  (2.15)  and  their 
boundary  conditions  can  be  written  in  the  form: 

•  1  de  ,  .  3F’  3f. 

(bp>)  *  5  *  x.  -  -  xir  f-  s> 

Y  3  0,  F  3  F'  =  0 


(2.36) 


Y  3  Ye,  F'  3  e,  e  -  cu  (Yee  -  W)  3  q± 


(2.37) 


(2.38a) 


(2.38b) 


whe  r  e 


Here  the  parameter  g^,  which  results  from  the  discrete  approximation  to  Eq. 
(2.19)  is  given  by 


where 


^  cij(Dj  '  Dj>  -  ctiDJ 


D  =  /vx/u  (Y  e  -  F  ) 
o  e  e 


(2.39) 


For  quasi- three-dimensional  flows,  we  define 


x  =  x, 


dn  =  (— )  dy 


(2.40) 


Here  u  is  a  reference  velocity  and  s,  denotes  the  length  in  the  longitud- 
o  1 

inal  direction  measured  from  the  initial  line  x  =  x*.  we  again  introduce  the 


definition  of  stream  function  f(x,y) 

df 

uh2  sln0  =  3y  - 


vV2sine  -  - 


(2.41)) 


so  that  with  the  definition  of  eddy  viscosity,  the  quasi-three-dimensional 
boundary- layer  equations  given  by  Eq.  (2.13)  can  be  written  as  [10] 

(bf-)'  +  ra^f  f  *  +  m3[(f’)2  -  e2]  +  m^f'g’  -  wge)  +  m4[(g')2  -  (we)2] 


/c,  ill  3e  ar 
=  Vf  aiT  “  e  3x  ~  £  3x} 


(2.42) 


(bg-)'  +  ra^f g"  +  ra^g')2  -  w2]  +  ra^f’)2  -  e2)]  +  nyf'g'  - 


,,,  iSLl  e  .  df. 

=  m9(f  3x”  "  e  aiT  ‘  9  3x) 

where,  with  f*  =  u/uo,  g'  =  w/Uq,  we  =  we/uQ 


(2.43) 


ml  =  hLh2  sin0  dx  (^1  h2  3ln0)' 


ra^  =  K^s^  cot© 


mj  =  ~K  s  csc6, 


-K,.,s ,, 


m,  =  s,K„  cot6 


(2.44) 


The  boundary  conditions  become 


n  =  0. 


f  -  f*  »  0. 


g  =  g'  =  0 


n 


f'  =  e(x). 


g’  =  «eU) 


~o 

w 

e 


where 


e(x) 


o.l  ,Xb  d&  do 

u  +  —  J  —  - 

e  *  1  do  x  -  o 

x 


f 

e 


). 


R 


u  s 
o 


V 


1 


(2.45a) 

(2.45b) 


(2.46) 


2.2.5  Solution  Procedure 

The  numerical  solution  of  the  system  of  equations  given  in  the  previous 
section  is  obtained  with  Keller's  box  method  for  the  standard  and  interactive 
methods.  This  is  an  efficient,  second-order  finite-difference  method  extens¬ 
ively  used  by  Cebeci  and  his  associates  for  a  wide  range  of  flows,  as  discussed 
in  Bradshaw  et  al.  [16].  The  description  of  the  standard  method  is  given  in 
that  reference  as  well  as  in  Cebeci  and  Bradshaw  [29].  The  general  features 
of  the  inverse  method  which  makes  use  of  the  Mechul- function  formulation  are 
also  described  for  two-dimensional  flows  in  Bradshaw  et  al.  [16]  and  in  [10] 
for  quasi-three-diraensional  flows.  As  in  previous  two-dimensional  studies  the 
FLAKE  approximation  in  which  the  convective  term  3f'/3x  is  set  equal  to 
zero  in  the  recirculating  region  is  employed,  and  no  attempt  was  made  to 
improve  the  accuracy  of  the  solutions  resulting  from  this  approximation. 


As  in  the  solution  of  two-dimensional  flows  by  Keller’s  method,  we  write 
Eq.  (2.37)  as  a  first-order  systtm.  For  this  purpose  we  let 


F*  =  u 
u'  =  v 


and  write  Eq.  (2.37)  as 


(bv)'  + 


1 

2 


F V  +  xe 


de 

dx 


,  3u  3f  , 

x(u  ^  -  v  ai  > 


(2.47a) 

(2.47b) 


(2.47c) 
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Since  e  Is  a  function  of  x.  only,  we  write 

e*  =  0  (2.47d) 

The  boundary  conditions  for  the  system  given  by  Bqs.  (2.47)  now  can  be  written 
as 

Y  =  0,  P  =  0 ,  u  =  0  (2.48a) 

Y  =  Y  .  u  =  e.  e  -  c .  (Ye-F)  *  g .  (2.48b) 

e  11  e  1 

After  the  finite-difference  approximations  to  Bqs.  (2.47)  and  (2.48)  ire 
written,  the  resulting  nonlinear  algebraic  system  Is  linearized  by  Newton’s 
method  and  the  linear  system  is  then  solved  by  the  block  elimination  method. 

For  further  details,  see  [29]. 


SaCTIOM  III 

BBULTS  MID  DISCUSSIONS 


N.V 
S'  V 


The  intet active  boundary* layer  procedure  described  In  the  previous 
section  has  been  applied  to  a  number  of  test  cases  and  the  results  Mill  be 
presented  here. 

Section  3.1  presents  the  results  obtained  frost  the  application  of  the 
two-dimensional  strip  theory  approach  to  two  wing  alone  configurations  to 
verity  s.e  soundness  of  the  Interactive  boundary  layer  technique.  Section  3.2 
will  present  the  results  obtained  frost  the  two-dimensional  and  the  quasi 
three  dimensional  methods  for  three  wing- body  fighter  configurations. 

3.1  iina  Mam  flaotAaiCiilaift 

The  first  conf lgurat Ion  considered  Is  an  RAB  clean  wing  with  a  28*  sweep 
angle  for  which  exper latent a 1  data  has  been  obtained  by  Lovell  [30].  Figure  5 
shows  the  computed  lift  curve  up  to  18*  angle  of  attack.  It  can  be  seen  that 
there  is  good  agreement  up  to  about  12*.  beyond  which  the  discrepancy 
Increases.  This  disagreement  at  the  higher  angles  of  attack  is  due  to  the 
three  dimensional  nature  of  the  flow,  and  due  to  the  sensitivity  of  the  flow 
to  the  transition  location.  This  sensitivity  Is  Illustrated  by  Fig.  6  which 
shows  the  variation  of  the  flow  solution  to  the  transition  location  at  an 
angle  of  attack  of  17.5*.  Moving  the  transition  location  leads  to  a 
significant  change  in  the  lift  separation  location  and  displacement  thickness, 
particularly  near  the  wing  tip.  Figure  7  shows  the  results  for  a  *  18.5* 
with  the  computed  transition  location  As  can  be  seen  the  separation  location 
at  this  angle  of  attack  occurs  at  about  the  midchord  location.  The  size  of 
the  leading  and  trailing  edge  separation  regions  at  a  =  17.5*  can  be  seen 
from  Fig.  8  in  which  the  local  skin  friction  Is  plotted. 

Mear  the  wing  tip  It  can  be  seen  that  there  Is  a  small  leading  edge 
separation,  shown  by  the  region  of  negative  skin  friction,  and  that  the 
trailing  edge  separation  Is  occurring  at  about  65%  local  chord. 

Figures  9  11  illustrate  the  results  obtained  for  a  MAOA  0012  swept  wing 
for  which  data  was  obtained  by  Yip  and  Shubert  [31]  for  a  range  of  angles  of 

f  hr  p  1  rjij  r m  Q  shCWS  *  y  \  1  nw  *  «;i  i  /i  lift  up  fn  71  f  f  \  <  h 

the  effects  of  the  Increasing  flow  separat Ion  at  higher  angles  of  attack  can 
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clearly  be  seen.  Figure  10  shows  the  computed  and  experimental  pressure  dis¬ 
tributions  for  a  -  19.35°  at  two  sections,  50%  and  85%  of  semispan,  both  of 
which  agree  very  well. 

The  Importance  of  the  transition  location  Is  Illustrated  by  Figure  11 
which  shows  the  computed  separation  location,  tralling-edge  displacement 
thickness  and  lift  distribution  for  a  =  21.12°.  Two  curves  are  shown,  the 
solid  line  was  obtained  with  a  computed  transition  location  at  about  3%  chord 
on  the  upper  surface.  The  dashed  line  was  obtained  by  specifying  a  transition 
location  closer  to  the  leading  edge  attachment  line  on  the  lower  surface.  It 
can  be  seen  that  this  small  movement  In  the  transition  location  moves  the 
separation  location  forward  from  60%  chord  to  about  30%  chord. 

3-2  application  to  Wing/ Body  Teat  Cases 

Three  wlng/body  test  cases  were  Identified  for  which  there  was 
experimental  data  available  to  evaluate  the  current  procedure.  The  three 
cases  are,  an  F-15  with  a  modified  wing  designed  to  Investigate  the  role  of  a 
leading- edge  laminar  separation,  an  Advanced  Navy  Fighter  configuration,  and 
an  unmodified  F15.  The  geometry  and  Input  data  associated  with  these  test 
cases  Is  discussed  In  detail  in  Section  5  (Figs.  26-30). 

The  F  15  laminar  bubble  configuration  has  been  run  at  a  range  of  angles 
of  attack  using  both  the  2D  and  the  quasi  3-D  boundary  layer  procedures. 
Figure  12  shows  the  computed  lift  compared  with  the  experimental  data  [32]. 

It  can  be  seen  that  at  lower  angles  of  attack  the  2D  and  quasi  3D  approaches 
agree  closely,  while  at  higher  angles  of  attack  the  quasi  3  D  boundary  layer 
procedure  agrees  better  with  the  experimental  data.  At  lower  angles  of  attack 
both  approaches  overpredlct  the  lift,  possibly  due  to  an  Inadequate  modelling 
of  the  flow  over  the  body. 

Figure  13  shows  an  example  of  the  computed  and  experimental  pressure 
distributions  for  a  =  10.75°  angle  of  attack  at  three  stations  across  the 
span  from  close  to  the  wing  tip  to  close  to  the  wing  root.  It  can  be  see., 
that  for  this  angle  of  attack  there  Is  very  little  difference  between  the  2  D 
and  the  quasi  three  dimensional  boundary  layer  methods,  and  that  both  give  a 
reasonable  approximation  to  the  measured  pressures  except  close  to  the  leading 
edge . 
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One  further  aspect  of  the  present  method  which  has  been  investigated 
concerns  the  use  of  a  panel  method  for  the  calculation  of  compressible  flow. 

A  panel  method  is  an  incompressible  solution  procedure  and  so  Mach  number 
effects  must  be  accounted  for  by  means  of  a  compressibility  correction.  In 
this  case  a  Goethert  correction  is  applied  using  the  following  procedure:  the 


solution  to  the  linearized  potential  flow  equation 

u  -  <>  &  +  .  &  .  0 

"  3x2  3y2  3z2 


(3.1) 


is  obtained,  where  $  is  the  perturbation  in  the  potential  due  to  the  body. 
Provided  that  this  perturbation  velocity  is  small  compared  with  the  freestream 
velocity  equation  (3.1)  provides  a  good  approximation  to  the  flow.  This  is 
true  for  the  flow  over  thin  wings.  However,  around  the  leading  edge,  or  on 
any  forward  facing  surface,  this  basic  assumption  is  violated. 

The  compressibility  correction  for  the  initial  velocity  calculation  is 

2 

Implemented  here  by  first  scaling  the  y  and  z  coordinates  by  3  »  V(1  -  M^) . 

The  Incompressible  flow  is  then  computed  about  this  equivalent  body  and  the 

2 

computed  velocity  components  are  divided  by  3  in  the  x-directlon  and  3 
in  the  y-  and  z-directions. 


For  the  subsequent  updates  to  the  inviscid  velocity  to  account  for  the 
viscous  effects  the  same  procedure  is  adopted.  However,  this  calculation 
involves  two  additional  features  over  the  initial  inviscid  calculation.  The 
first  is  the  introduction  of  the  blowing  velocity  required  to  simulate  the 
displacement  thickness  and  the  second  is  the  introduction  of  off-body  points 
at  which  the  Kutta  condition  is  to  be  applied  and  the  final  velocity  computed. 
In  applying  the  Mach  number  correction  two  additional  assumptions  are  there¬ 
fore  necessary.  The  first  is  that  the  blowing  velocity  is  predominantly  in  the 
y.z  direction  so  that  the  blowing  velocity  applied  on  the  "equivalent  incom¬ 
pressible  body"  can  be  scaled  by  3.  The  second  assumption  is  that  the  dis¬ 
placement  effect  of  the  boundary  layer  is  also  predominantly  in  the  y,z  direc¬ 
tion.  Therefore  in  defining  the  off-body  point  location  in  the  transformed 
plane  the  displacement  thickness  is  first  scaled  by  3.  Since  both  the 
displacement  thickness  and  the  blowing  velocity  are  most  significant  near  the 
trailing  edge  both  of  these  assumptions  are  valid. 


It  should  be  pointed  out  that  the  Interactive  viscous  calculation  is 
performed  entirely  In  the  physical  plane.  The  Mach  number  correction  outlined 
above  Is  used  to  provide  the  compressible  external  velocity  distribution. 

While  the  compressibility  correction  will  provide  a  good  approximation  at 
lower  angles  of  attack,  it  becomes  less  accurate  at  higher  angles  of  attack 
when  the  flow  can  become  supercritical  and  shock  waves  can  form  through  the 
leading  edge  region.  Figure  14  shows  the  inviscid  computed  pressure 
distribution  at  the  mid-semispan  location  computed  by  the  present  panel  method 
and  by  the  transonic  full-potential  wing/body  code  developed  by  Chen  et  al. 

[33]  at  a  =  6.84°,  10.75°  and  12.95°  for  M  =  0.6.  At  the  lower  angle  of 
attack  the  two  procedures  agree  well  except  close  to  the  leading  edge  where 
the  transonic  code  predicts  a  higher  suction  peak.  However,  at  a  =  12.95° 
the  transonic  method  predicts  a  shock  wave  and  a  larger  difference  in  the 
pressure  distribution  through  this  region  is  predicted. 

The  second  wlng/body  test  case  which  has  been  considered  is  the  Advanced 
Navy  Fighter  (ANF)  for  which  experimental  data  is  available  [34].  This  con¬ 
figuration  was  tested  experimentally  both  with  and  without  a  canard.  Figure  15 
shows  the  computed  and  experimental  lift  distribution  for  this  configuration 
without  the  canard  and  Figure  16  shows  the  viscous  and  inviscid  pressure  dis¬ 
tribution  compared  with  experimental  data  at  o  =  7.7°  across  the  span.  It 
can  be  seen  that  there  is  good  agreement.  The  only  other  pressure  data  for 
this  particular  configuration  is  for  a  =  22°  at  which  angle  the  data  shows 
the  presence  of  a  shock  wave  near  the  leading  edge  suction  peak.  Due  to  the 
limitations  of  the  compressibility  correction  outlined  above  this  case  has  not 
been  pursued.  The  effect  of  the  canard  is  shown  in  Figure  17,  in  which  the 
computed  viscous  lift  distribution  is  shown  with  and  without  the  canard.  This 
case  was  run  with  the  canard  at  zero  deflection  angle  for  which  there  was  no 
experimental  data,  and  the  viscous  effects  were  computed  only  on  the  main  wing. 

The  final  configuration  which  has  been  investigated  is  an  unmodified  F- 15 
wing/body  for  which  experimental  data  was  obtained  by  Anderson  [35].  Figure  18 
shows  the  computed  and  experimental  lift  curves  from  which  it  can  be  seen  that 
there  is  good  agreement  up  to  8.66°  angle  of  attack.  Figure  19  shows  the 
computed  pressure  distribution  for  8.66°  which  is  in  good  agreement  with  the 
experimental  data. 
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SECTION  IV 
PROGRAM  DESCRIPTION 


The  method  developed  here  calculates  the  viscous  flow  over  a  lifting  body 
through  the  use  of  an  interactive  boundary  layer  technique  In  which  a  three 
dimensional  panel  method  Is  coupled  with  a  finite-difference  boundary  layer 
method.  Two  distinct  boundary- layer  procedures  have  been  adopted,  the  first 
of  which  uses  a  two-dimensional  strip  theory  approach  while  the  second  Is  a 
quasi-three- dimensional  boundary- layer  method. 

The  structure  of  the  code  is  illustrated  in  Figure  20,  which  gives  a 
schematic  flow  chart  of  the  code.  There  are  three  major  computational  modules 
required,  namely  the  inviscld  panel  method,  the  two-dimensional  boundary  layer 
routines  and  the  quasi- three-dimensional  boundary- layer  routines.  The  commun 
lcation  between  the  potential  flow  code  and  the  boundary- layer  routines  Is 
handled  through  the  use  of  external  disk  storage. 

For  a  given  flow  condition  the  inviscld  solution  Is  first  computed  and  the 
pressure  distribution  over  the  wing  is  saved.  If  the  2d  boundary- layer  mode 
is  selected,  this  data  is  separated  into  strips,  and  an  interactive  boundary- 
layer  calculation  is  performed  independently  for  the  upper  and  lower  parts  of 
each  lifting  strip  using  the  method  outlined  in  Section  2.2.5.  Each  of  these 
boundary- layer  calculations  require  a  number  of  iterative  sweeps  over  each 
surface  In  order  to  match  the  current  external  velocity  distribution,  given  by 
Eqs.  (2.18)  and  (2.19),  with  that  computed  by  the  boundary  layer  equations, 

(2.47)  and  (2.48).  Once  these  inner  Iterations  have  been  completed  and  the 
corresponding  displacement  thickness  and  blowing  velocity  has  been  computed 
these  data  are  transferred  back  to  the  Inviscld  flow  code  and  the  external 
velocity  distribution  Is  updated  to  incorporate  the  viscous  effects.  One 

complete  Iteration  between  the  viscous  and  Inviscld  codes  Is  referred  to  as  a  < 

\ 

cycle.  Several  such  cycles  are  usually  performed  in  order  to  obtain  a  fully  j 

converged  solution  In  which  the  computed  viscous  and  Inviscld  solutions  match  I 

one  another.  The  precise  number  of  cycles  performed  is  governed  either  by  an  ' 

Input  parameter  which  specifies  the  maximum  number  of  cycles,  or  else  by  a 
convergence  check  which  Is  based  on  the  maximum  change  In  a  displacement  j 

thickness  parameter  form  one  cycle  to  the  next.  Typically  up  to  5  cycles  may  I 

be  required  for  a  fully  converged  solution.  ! 
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A  similar  procedure  employed  for  the  quasi  three  dimensional  boundary 
layer  method  with  interface  routines  handling  the  transfer  of  the  data  between 
the  viscous  and  the  invlscid  codes.  In  this  case  the  data  required  is  more 
extensive  than  that  required  for  the  two-dimensional  boundary  layer  calculation 
since  surface  coordinates  and  curvatures  are  required  for  the  complete  surface. 
Therefore,  the  interface  program  accesses  the  velocity  data  produced  by  the 
invlscid  calculations,  and  computes  a  dividing  line  separating  the  upper  and 
lower  portions  of  the  wing  flowfield.  For  each  surface,  a  boundary  layer  grid 
Is  now  defined  and  the  geometric  surface  curvatures  and  the  velocity  components 
ate  Interpolated  onto  the  boundary  layer  grid.  The  quasi  three  dimensional 
boundary  layer  procedure  now  Iterates  over  each  surface  in  a  similar  fashion 
to  the  two  dimensional  approach  in  which  the  perturbed  invlscid  velocity  is 
matched  with  that  computed  by  the  boundary  layer  equations,  (2.42)  (2.46). 

The  computed  displacement  thickness  and  blowing  velocity  is  then  transferred 
back  to  the  potential  flow  code  for  the  next  viscous/ invlscid  cycle. 

After  a  converged  solution  has  been  obtained  for  the  first  specified  angle 
of  attack,  the  process  is  repeated  for  the  next  angle  of  attack.  However, 
rather  than  start  the  calculation  using  the  purely  invlscid  potential  flow 
solution,  the  blowing  velocity  and  displacement  thickness  computed  for  the 
previous  angle  of  attack  are  used  for  boundary  conditions  in  the  first  invlscid 
flow  solution.  Therefore,  when  running  multiple  angles  of  attack,  the  angles 
should  be  specified  in  increasing  order.  In  this  way  the  number  of  cycles 
required  for  a  converged  solution  is  reduced,  particularly  at  higher  angles  of 
attack  when  boundary  layer  separation  can  lead  to  fairly  large  displacement 
thicknesses.  To  take  further  advantage  of  this  fact,  the  program  may  be  exe 
cuted  In  a  restart  mode  in  which  data  stored  by  a  previous  run  may  be  input  as 
a  starting  solution  for  the  next  angle  of  attack  to  be  computed.  After  each 
run,  the  following  two  datasets  may  therefore  be  saved  for  use  in  a  subsequent 
run : 


Unit  #33 


and 

Unit  #47 


This  stores  the  blowing  velocity  and  the  displacement 
thickness  data. 

This  unit  stores  the  geometric  coefficients  of  the  airfoil 
appearing  in  the  interaction  formula  used  In  the  quasi 
three  dimensional  boundary  layer  method. 


2b 
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Unit  #51  -  This  unit  stores  the  geometric  coefficients  of  the  air¬ 
foil  appearing  in  the  interaction  formula  used  in  the 
two-dimensional  boundary- layer  method. 

The  user  can  therefore  run  a  range  of  angles  of  attack  either  by  specify¬ 
ing  all  of  the  angles  of  attack  in  one  single  run,  or  else  as  a  series  of 
separate  runs  which  make  use  of  the  restart  capability. 

4.1  Input  P«u  Description 

4.1.1  Body  Oeoaatry 

The  input  to  this  program  consists  of  the  coordinates  of  a  number  of 
points  that  define  the  surface  of  a  three-dimensional  configuration  on  which 
the  flow  is  to  be  computed.  For  the  purpose  of  organizing  these  points  for 
computation,  each  point  is  assigned  a  pair  of  integers,  m  and  n.  These 
integers  need  not  be  input,  but  their  use  must  be  understood  to  insure  the 
correctness  of  the  input  and  to  facilitate  the  interpretation  of  the  output. 

For  each  point,  n  identifies  the  "column*  of  points  to  which  it  belongs, 
while  m  identifies  its  position  in  the  "column,"  i.e,  the  "row."  The  first 

point  of  a  "column*  always  has  m  ■  1.  To  insure  that  the  program  will  compute 

outward  normal  vectors,  the  following  condition  must  be  satisfied  by  the  input 
points.  If  an  observer  is  located  in  the  flow  and  is  oriented  so  that  locally 

he  sees  points  on  the  surface  with  m  values  increasing  upward,  he  must  also  see 

n  values  increasing  toward  the  right.  Examples  of  correct  and  incorrect  input 
are  shown  in  Plgure  21.  In  this  figure  the  flowfield  lies  above  the  paper, 
while  the  interior  of  the  body  lies  below  the  paper.  Occasionally,  it  happens 
that  despite  all  care  a  body  is  input  incorrectly.  If  the  entire  body  is  input 
incorrectly  -  not  some  sections  correctly  and  some  incorrectly  -  the  difficulty 
can  be  remedied  by  changing  the  sign  of  one  coordinate  of  all  the  input  points. 
This  trick  will  give  an  input  body  of  the  proper  shape  at  perhaps  a  peculiar 
location.  Otherwise,  the  input  will  have  to  be  done  over. 

The  body  surface  is  divided  into  sections,  which  may  be  actual  physical 
divisions  or  may  be  selected  for  convenience.  A  section  is  defined  as  con¬ 
sisting  of  a  group  of  at  least  two  n-lines.  Within  each  section  the  n-llnes 
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are  input  in  order  of  increasing  n.  On  each  n-line  the  points  are  input  in 
order  of  increasing  ra.  All  n-lines  in  a  section  must  have  the  same  number  of 
points,  but  this  may  vary  from  section  to  section.  The  first  n-line  of  the 
first  section  is  n  =  1.  Prom  then  on  the  n-lines  may  be  thought  of  as  numbered 
consecutively  through  all  sections,  i.e.,  the  numbering  is  not  begun  over  at 
the  beginning  of  each  section.  Panels  will  be  formed  that  are  associated  with 
points  on  every  n-line  except  those  that  are  last  in  their  respective  sections. 
Points  on  these  latter  n-lines  are  used  only  to  form  panels  associated  with 
points  on  the  next  lowest  n-lines. 

To  illustrate  this  procedure,  consider  the  plan  view  of  a  body  shown  in 
Figure  22.  This  body  has  been  divided  into  four  sections,  as  shown  in  the 
figure.  The  first  section  contains  tour  n-lines,  n  =  1,  2,  3,  4;  the  second, 
five  n-lines,  n  =  5,  6,  7,  8,  9;  the  third  three  n-lines,  n  =  10,  11,  12;  and 
the  fourth  three  n-lines;  n  =  13,  14,  15.  The  number  of  points  on  each  n-line 
are: 


Section  =1234 
M  =  4  7  4  2 

Notice  that  the  line  n  =  4  has  only  four  points,  the  points  m  =  1,  2,  3,  4  and 
the  ra-grid  of  Section  1,  which  is  listed  in  the  figure  along  the  n  =  1  line. 

The  lines  n  =  4  and  n  =  5  are  physically  identical,  some  of  the  points  on  the 
two  lines  are  physically  identical  but  correspond  to  different  values  of  m. 

This  is  of  no  consequence.  In  this  scheme  sections  are  completely  independent. 
No  elements  are  computed  corresponding  to  points  on  lines  n  =  4,  9,  12,  15. 

There  is  no  restriction  tha:  the  m  and  n  lines  of  different  sections  have 
to  be  roughly  parallel.  The  arrangement  shown  in  Figure  23  is  permissible. 

As  discussed  in  Section  2.1,  the  body  is  divided  into  lifting  and  nonlift¬ 
ing  sections.  The  arrangement  of  the  input  requires  that  all  lifting  sections 
precede  nonlifting  sections  with  a  one-point-per-card  format  of  (3F10.5,2I1) 
with  a  maximum  of  2000  panels  and  100  strips  defining  the  configuration.  This 
data  is  input  to  the  program  through  Unit  #29.  Also,  since  the  boundary- layer 
calculations  can  be  made  for  only  one  lifting  section,  this  section  must 
precede  any  other  lifting  sections.  The  two  integers  that  follow  each  set  of 
coordinates  are  the  status  flag  and  label  flag,  respectively,  that  is 
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STATUS  FLAG  LABEL  FLAG 

2  -  new  section  1  -  lifting  body,  or  0  -  nonllftlng  body 

l  -  new  N-line  (these  flags  are  needed  only  when  STATUS  FLAG  =  2) 

0  -  same  N-line 

3  -  end  of  input 

for  the  whole  body 

4.1.2  Namelist  Block  Input 

There  are  two  namelist  blocks  that  contain  two  groups  of  Input  data. 


Namelist  Name  A: 


This  namelist  block  consists  of  all  the  control  flags  and  flow  conditions 
required  to  control  the  interactive  viscous  calculation. 


Variables  Remarks 

LIFT3D  Control  flag  for  the  execution  of  the  potential-flow 

program. 

=  0  Skip  the  Initial  potential  flow  and  execution  will 
start  with  the  boundary- layer  code. 

=  1  Use  the  potential  flow  program  for  every  cycle. 

(default  is  1) 

IBL3D  Control  flag  for  executing  a  boundary-layer  program. 

=  0  Use  the  two-dimensional  boundary- layer  program. 

=  1  Use  the  quasi-three-dimensional  boundary- layer  program, 
(default  is  0) 

KINKNO  Total  number  of  "kink"  stations  input  ("kink"  station  means 

a  discontinuity  in  the  planform  of  the  first  lifting 
section) . 

(default  is  0,  maximum  is  5) 

MS  Station  numbers  where  "kink*  occurs  in  the  planform.  The 

order  of  these  "kink"  stations  should  be  in  the  same  order 
as  the  spanwise  stations  input  to  the  Neumann  program. 

They  are  also  part  of  the  input  spanwise  stations. 

(default  is  0) 

MLINE(i)  Number  of  defining  chordwise  points  on  each  N-line  of  the 

i-th  lifting  section. 

NLINE(l)  Number  of  defining  spanwise  points  on  each  M-line  of  the 

i-th  lifting  section. 

MAXCYC  Maximum  number  of  interaction  cycles. 

(default  is  0) 
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Variables 


Remarks 


t 

s 


>r* 

Sf 

S 


% 

%■ 

V 


v 
■  ^ 

V 

v 


MAXANG 

IPNTGM 

IPNTNU 

IPNTB3 

DBMAX 

ALPHA I 

AMACH 


Number  of  angles  of  attack  to  be  calculated. 

(default  is  1,  maximum  is  10) 

Optional  (debug)  print  flag  for  the  interfacing  program. 

=  0  No  intermediate  print. 

=  1  With  intermediate  print. 

(default  is  0) 

Optional  (debug)  print  flag  for  the  potential  flow  program. 
=  0  No  intermediate  print. 

=  1  With  intermediate  print. 

(default  is  0) 

Optional  (debug)  print  flag  for  the  quasi-three-dimensional 
boundary  layer  program. 

=  0  No  intermediate  print. 

=  1  With  intermediate  print. 

(default  is  0) 

Maximum  difference  in  circulation  values  generated  from 
potential  flow  solution  for  convergence. 

(default  is  0.01) 

Input  angles  of  attack  values. 

Total  number  input  =  MAXANG. 

(default  is  0.0) 

Input  Mach  number. 

(default  is  0.0) 


Namelist  Name  B: 

Namelist  B  contains  the  geometric  variables  and  control  flags  which 
pertain  to  the  potential  flow  calculation. 


Variables  Remarks 

BOV 2  The  semi-span  value  of  the  lifting  body, 

(default  is  1.0) 


IAUTOW  Automatic  wake  generation  flag: 

=  1  Program  generates  bisector  wake. 

=  2  Program  generates  flat  wake,  i.e.  wake  panel  is 
parallel  to  the  x-axis. 

=  0  User  input  wake  panels. 

(default  is  1) 

IPCV  Chordwise  vorticity  flag. 

=  0  Constant  vorticity  used. 

=  1  Parabolic  vorticity  used. 

(default  is  0) 


> 
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Variab les 


Remarks 


ITERAT 


IVNSP 


LIFSEC 


LINEAR 


NOTE:  This  flag  controls  the  form  of  the  predetermined 
chordwlse  bound-vorticity  variation  discussed  in  Section 
2.1.  For  wings  with  cuspid  or  thin  trailing  edges,  the 
constant  variation  may  lead  to  unrealistic  variations  in 
the  pressure  close  to  the  trailing  edge.  In  such  cases, 
the  parabolic  variation  should  be  used.  Otherwise  the 
constant  chordwlse  option  should  be  used. 

Matrix  Solution  flag. 

=  0  (direct  solution  for  350  panels  or  less,  otherwise 
use  iterative  solution). 

=  1  Iterative  solution. 

=  2  Direct  solution. 

(default  is  1) 

Note:  Iterative  solution  must  be  used  when  viscous 
calculations  are  made. 

Nonuniforra  Flow/Specified  Normal  Velocity  flag. 

=  0  No  nonuniform  flow  or  blowing  velocity. 

=  1  Nonuniform  Flow  specified.  In  the  case  onset  flow 
velocity  components  at  every  control  point  must  be 
input  on  Unit  33  in  3F10.6  card  image  format. 

Otherwise  velocity  components  are  computed  from 
freestream  flow  direction. 

=  2  Nonzero  normal  velocity  condition  is  specified  on  some 
section  or  sections.  If  this  option  is  used,  then 
Unit  33  must  contain  the  following  data: 

1st  record:  "SECTION  k" 

Subsequent  records:  required  normal  velocity  (FORMAT 
1F10.6)  with  1  record  for  each  control  point  in  the 
k-th  section.  These  cards  are  repeated  for  each 
section  for  which  nonzero  boundary  condition  is 
satisfied. 

NOTE:  These  options  are  only  available  for  a  purely 
inviscid  calculation. 

Number  of  lifting  sections  input. 

NOTE:  Boundary-layer  calculations  can  be  made  for  the 
first  lifting  body  only. 

Spanwise  vorticity  Elag. 

=  0  Constant  vorticity  used. 

=  1  Linear  vorticity  used.  If  selected,  the  variables, 
NLINE1  and  NLINEN  must  be  input. 

(default  is  0) 

NOTE:  This  flag  controls  the  variation  of  spanwise 
vorticity  across  a  given  lifting  strip.  The  level  on  each 
strip  is  determined  by  the  Kutta  condition  as  discussed  in 
Section  2.1. 


Variables 

NLINEl(i) 


NLINEN( i) 


LIST 


Remarks 


Condition  for  the  first  N-line  of  the  i-th  lifting  section. 
=  1  The  first  N-line  is  the  wing  tip. 

=  4  The  first  N-line  belongs  to  an  extra  strip  or  is  next 

to  a  symmetry  plane. 

Condition  for  the  last  >l-line  of  the  i-th  lifting  section. 

=  1  The  last  N-line  is  the  wing  tip. 

=  4  The  last  N-line  belongs  to  an  extra  strip  or  next  to  a 

symmetry  plane. 

Geometry  formation  flag. 

=  0  Complete  potential  flow  calculation  modes. 

=  1  Calculation  stop  after  geometry  formation 
(default  is  0) 


NOFF  Off-body  point  flag. 

=0  No  off-body  point  input. 

=  1  Off-body  point  input. 

(default  is  0) 

NOTE:  Off-body  points  are  points  in  the  flow  at  which  the 
velocity  and  pressure  are  to  be  computed.  If  this  option 
is  used,  the  off-body  coordinates  must  be  supplied  on  Unit 
32  using  3F10.6  card  image  format. 

NSYM1  First  plane  of  symmetry  flag. 

=  0  No  symmetry. 

=  1  One  plane  of  symmetry  (about  the  x-z  plane). 

(default  is  1) 


NSYM2  Second  plane  of  symmetry  flag. 

=0  No  second  plane  of  symmetry  used. 

=  1  Symmetry  about  the  y-z  plane  (NSYMl  must  also  be  =  1). 
(default  is  0) 


v-v 

A" 

ORIGNX 

ORIGNY 

ORIGNZ 

NSORCE(i) 

m 

VN 

fXv- 

NWAKE ( i ) 

W 

,  r 

NSTRIP(i) 

,v.v 

IXFLAG(i) 

•>  V 

pi 
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The  x.y.z  coordinates  of  the  input  moment  origin 
(default  is  0,  0,  0). 


Number  of  on-body  panels  per  strip  for  the  i-th  lifting 
section. 


Total  number  of  wake  panels  input  for  the  i-th  lifting 
section. 

(If  IAUTOW  *  0,  NWAKE  5  1) 


Total  number  of  strips  in  the  i-th  lifting  section. 


Extra  strip  option  input  for  the  i-th  lifting  section. 
=  0  No  extra  strip. 

=  i  First  strip  is  an  extra  strip. 

=  2  First  and  last  strips  are  extra  strips. 

=  3  Last  strip  is  an  extra  strip. 


32 


Variables 


Remarks 


NOTE:  An  extra  strip  is  a  strip  of  panels  on  a  lifting 
section  which  carries  a  vorticity  distribution  but  on  which 
no  boundary  conditions  are  satisfied.  Such  a  strip  is  used 
for  instance  at  an  intersection  between  a  wing  and  a 
fuselage  in  order  to  provide  a  more  accurate  carry-over  of 
the  wing  lift  across  the  fuselage.  In  this  case  the  extra 
strip  is  internal  to  the  fuselage,  and  it  extends  from  the 
wing  root  to  the  fuselage  centerline.  Further  details  of 
the  use  of  extra  strips  is  discussed  by  Hess  [36]. 

REFCHD  Reference  chord  value.  (This  is  needed  for  the  calculation 

of  the  lift  coefficient.) 

(default  is  1.0) 

RFAREA  Reference  area  value.  (This  is  needed  for  the  calculation 

of  the  lift  coefficient.) 

(default  is  1.0) 


4.1.3  Boundary-Layer  Calculation  Data  Input 

(a)  Card  number  1 
Format  (6F10.0) 

Card  Column  Variables  Remarks 

1-10  RCxlO"6  Reynolds  number  based  on  CREF. 

11-20  CREF  Reference  length.  In  the  same  physical  units 

that  are  used  for  the  input  geometry. 

NOTE:  This  need  not  be  the  same  value  as  REFCHD  used 
above.  CREF  is  only  used  to  calculate  the  local 
Reynolds  number  in  terms  of  the  input  coordinates, 
e.g.  If  the  geometry  is  specified  in  inches,  and  the 
Reynolds  number/ft  is  input,  then  CREF  would  be  12. 


31-40  FKSTAT  =  1.0  Option  for  using  a  previous  angle  of 

attack  solution  as  an  initial 
estimate  for  the  current  angle  of 
attack  calculation. 

=  0.0  Option  not  used  (or  the  initial  run). 

41-50  FPRNT  =  1.0  Detailed  boundary- layer  print. 

=  0.0  No  boundary- layer  print. 
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(b)  Card  Number  2 
Format  (6P10.0) 

Card  Column  Variables  Remarks 

1-10  FNBLU  Number  of  boundary  layer  stations  at  which 

calculations  are  made  on  the  upper  surface  of 
the  first  lifting  body  (max  -*  100). 

11-20  FTRLU  Transition  Input  flag 

-  0  Transition  locat Ion  specified  across 
the  span.  See  Cards  3  and  4. 

=  1  Transition  calculated  by  the  program  at 
pressure  peak. 

=  2  Transition  calculated  by  the  program 
using  Michel's  criteria. 

21-30  FSVPU  Number  of  Iteration  sweeps  to  be  made  for  the 

upper  surface. 

( FSVPU  T  10.0  is  recommended ) 

31~40  FNTRU  Nixober  of  transition  locations  to  be 

specif  led . 

(if  FTRIU  =  0) 


(c)  Card  number  3  (needed  only  if  FTRIU  =0). 

Format  (6F10.0) 

Card  Column  Variables  Remarks 

1-10  YTRN  Spanwlse  location.  In  the  input  geometry 

system,  where  transition  will  be  specified 
(this  location  does  not  have  to  correspond  to 
the  coordinates  used  to  define  the  geometry). 

YTRN ♦  1 


YTRN ♦ 5 

(Repeat  this  card  If  necessary  until  l  -  FNTRU.) 


(d)  Card  number  4  (needed  only  If  FTRlue  0). 

Format  (6F10.0) 

Card  Column  Variables  Remarks 

1_10  XCTRN  x/c  chordwlse  location  of  spe<  If  led 

t  r  ans  It  1  >n . 


11-20 


51-60 
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XCTRM+ 1 


5 1 -60  XCTRM+5 

(Repeat  this  card  If  necessary  until  1  *  FHTRU.) 

The  following  cards,  numbers  5,  6  and  7,  have  the  same  format  as  cards 
2.  3  and  4  and  the  variables  have  the  same  meaning  but  correspond  to  the  lower 

surface . 

le)  Card  number  5 
Format  (6F10.0) 

variables  FNBLL,  FTRIL,  FSWPL,  PNTRL 
FSWPL  =  5.0  Is  recommended. 

( f )  Card  niaber  6 
Format  (6F10.0) 

Variable  YTRN 

(g)  Card  number  7 
Format  (6F10.0) 

Variable  XCTRN 

Note  that  the  chordwlse  location  of  specified  transition  must  be  preceded 
by  a  minus  sign  when  It  Is  being  Input  on  the  airfoil  lower  surface. 

The  cards  In  this  section,  (a)  through  (g),  are  to  be  repeated  for  each 
angle  of  attack  calculated. 


<  2  gA&gai-Btu  BM£ilBiian 

4.2.1  Summary 
Potential  Flow 

The  following  data  are  printed  for  the  potential  flow  calculation: 
i .  Input  f lags . 


15 
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Tl.  Panel  formation. 


III.  Matrix  solution  (by  Iteration  method). 


IV.  Solution  scmnary.  This  Includes: 

lift  coefficient,  drag  coefficient,  pitch,  roll,  yaw  and  n 
values  for  each  lifting  and  nonllfting  section,  plus  the  same 
set  of  data  for  the  whole  configuration. 


v.  Final  cycle.  For  every  panel,  the  following  data  are  printed: 

.  panel  number 

.  control  point  coordinates  -  x  ,  y  ,  z 

o  o  o 

.  lnvlscld  velocities  -  v  ,  v  ,  v 

x  y  z 

.  lnvlscld  total  velocity  and  viscous  corrected  pressure 

coefficient  Cp 

.  displacement  thickness  DELS 
.  shape  factor  -  H 
.  skin  friction  -  CF. 


4.2.2  Format 

A  test  case  consisting  of  one  lifting  and  four  nonllfting  sections  with  a 
total  of  221  panels,  Figure  24,  was  created  to  present  a  sample  of  the  program 
output,  and  a  listing  of  this  geometry  (Input  through  Unit  #29)  Is  presented 
In  Table  1.  Two  angles  of  attack  were  calculated  with  the  transition  locations 
on  the  lifting  section  computed  by  the  program.  The  namelist  block  Input  for 
this  case  Is: 

REPORT  TEST  CASE, TRANS  CALC.,  2DSTRIP  BL,  H=0  8,  ALPHAS 
AA  LIFT3D=1, IBL3D=0, 

0AXCYC=2,  NLINE=5,  f1LINE  =  21 , 

0AXANG=2, ALPHAI (1) si . 3, ALPHAI (2):3  4, 

AnACHsO.e, AEND 

«B  B0Y2  =  12 . 955,  I  All  TOW:  2,  NSYfll  =  1 ,  LIFSECsi,  REFCHD=  8  2912 
UNEARsO,  IXFLAG(l)  s3#  JPCVsl, 

RFAREAsBS . 904,  NS0RCE=20, NSTRIP=4,  NWAKE  =  1, AEND 


7.0 

12.0 

0.0 

0.0 

70.0 

2.0 

10.0 

51.0 

2.0 

5.0 

7.0 

12.0 

0.0 

0.0 

70.0 

2.0 

10.0 

51 .0 

2.0 

5.0 

For  conciseness,  the  output  presented  here  will  be  the  pertinent 
information  for  the  first  angle  of  attack  only  and  Is  shown  in  Figure  2b. 
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2*. 320007  12. 260000  6. ••9997221  22.049991  6.9999999  6.479327200  26.103776  0.0000000  9 .995799300 
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4.3 


There  are  a  total  of  50  external  units  required  although  they  are  not 
simultaneously  used.  A  list  of  the  unit  numbers  and  types  follows: 

1.  Direct  access  data  sets: 

Unit  #20  -  400  records,  each  200  words  long  -  potential  flow 
velocity 

Unit  #61  -  3280  records,  each  990  words  long  -  used  for  quasi- 
3-D  boundary- layer  program 

Unit  #62  -  41  records,  each  990  words  long  -  used  for  quasi- 3-D 
boundary- layer  program 

2.  Sequential  data  sets:  This  Includes  units  3,  4,  5,  8,  9,  10,  11, 

12,  13,  14.  15,  16,  17,  18,  19,  22,  23,  24,  25,  26,  27,  28,  29,  30, 

31,  32,  33,  34,  35,  36,  38,  39,  40.  41,  42.  43.  44.  45,  56.  47.  51. 

63,  64.  65,  66,  67. 

Unit  #29  is  the  geometry  data  set  for  the  Input  body.  These 
datasets  are  used  both  for  the  lnvlscid  velocity  calculation  and  for 
the  communication  between  the  viscous  and  the  lnvlscid  parts  of  the 
calculation. 

3.  Datasets  required  to  be  saved: 

Unit  #34  -  This  unit  contains  the  lnvlscid  velocities  and 
control  points. 

Unit  #33  -  This  unit  contains  the  blowing  velocities  and  the 
displacement  thickness  values  for  viscous/ lnvlscid 
Interactive  runs. 

Unit  #51  -  This  unit  contains  some  input  data  to  the  2-D 
boundary- layer  program. 

Unit  #47  -  This  unit  contains  the  blowing  velocities  from  the 
3-D  boundary- layer  program. 

4.4  CRAY  JCL 

An  example  of  the  JCL  required  for  the  execution  of  the  program  on  CRAY 
using  the  two-dimensional  strip  theory  boundary- layer  procedure,  follows. 
Example  1  Is  for  an  Initial  calculation  with  the  appropriate  data  saved  for 
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subsequent  runs  and  example  2  shows  how  the  previously  saved  data  Is  used  to 
make  additional  calculations. 


Example  1 


JOB , JN  =  xxxxx  ,T=900,MFL=2000000  . 

ACCOUNT , US= xxxxxx  ,UPW  =  xxxxx  ,AC=XXXXXXXX 
« . 

COPYF , I=$IN,0=I NPUT .  COPY  INPUT  FILE 
REWIND, DN=INPUT. 

*  . 

COPYF, I=$IN,0=GE0M.  COPY  INPUT  FILE 
REWIND, DN=GE0M. 

* 

ACCESS, DN=KCMCAIR,PDN=KCMCAIRCRAY, ID=0BJLIB. 

STATIC, LEVEL=NEW. 

BUILD, 1=0, L=0,OBL=kCMCAIR, B=0 , NBL =MCARBI N , NODIR .  CREATE  SEQUENTIAL 

* 

ASSIGN, DN=INPUT,A=FT05.  ALLOCATE  APPROPRIATE  FILES 
ASSIGN, DN=$OUT,A=FT06  . 

ASSIGN, DN=GEOM, A =FT29. 

ASSIGN, DN=FILE34,A=FT34. 

ASSIGN, DN=FILE33,A=FT 33. 

ASSIGN, DN=FILE51,A=FT51 . 

*.  . 

MODE , FI =DISABLE . 

* . 

SEGL DR , CMD= ' B I N=MCARBIN ' .  LINK  BINARIES 
SABD  EXECUTE  LOAD  MODULE 


SAVE,DN=FILE34,PDN=INVELNEW.  SAVE  OUTPUT  DATASET  FOR  THE  NEXT  CAL. 
SAVE,DN=FILE33,PDN=BL  0W2DNEW . 

SAVE ,DN=FILE51,PDN=BL2DINPNEW. 

« 

EXIT  . 

/EOF 

FI  5  WITH  LAMINAR  BUBBLE  WING,  MULTIPLE  ALPHAS,  M=0.6 
BA  LIFT3D=1,IBL3D=0, MAXCYC=2 ,MLINE=51,NLINE=14, 

MAXANG=2, ALPHA  I ( 1 )=6 .95, ALPHAI(2)=8 .95, AMACHI 1 )=0 .60, 

AMACH(2)=0  6, SEND 

SB  BOV 2=13. 008, IAUTOW=2, NSYM1 =1,LIFSEC=1,NSTRIP=13, 

IXFLAGt 1 )=3,NSORCE=50,NWAKE=1 . RFAREA=96  . 7,REFCHD=8 . 325, SEND 


6  022 

8 . 325 

0 . 0 

0 . 0 

61.0 

2 . 0 

10.0 

0 . 0 

51  .  0 

2  0 

5.0 

0 . 0 

6  022 

8 . 325 

1  .  0 

0  .  0 

610 

2  0 

10  0 

0 . 0 

51  0 

2  0 

5.0 

0 . 0 

/EOF 

32  966003 
924698 
863266 
782425 
683380 
567581 


32 

32 

32 

32 

32 


13.008000 
13.008000 
13 . 008000 
13.008000 
13.008000 
13 . 008000 


949856821 

950995400 

952991500 

956209200 

959699600 

961723300 


43.000000  1.3186102 
43  000000  0 . 5522600 
43  000000  0  1862400 
43  000000  0.0000000 
/EOF 


5  600000400 
5 .600000400 
5  600000400 
5  600000400 


BIN 


saved  for 
-  future 
runs 


input 
data  set 


y  geometry 
data  set 
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JOB . JN=«xxxxx  ,T=1200,MFL= 2400000. 

ACCOUNT, US=xxxxxxx  ,UPW=XX**XXX,AC=XXXXXXXXX  . 

< . 

COPYF . I =$IN, 0= INPUT .  COPY  INPUT  FILE 
REWIND, DN=INPUT. 
x . 

COPYF , I =$IN, 0=GE0M.  COPY  INPUT  FILE 
REWIND, DN=GEOM. 

».  ACCESS  PERMANENT  DATASETS 
ACCESS, DN-KCMCAIR, PDN=KCMCAIRCRAY, ID=OBJLIB. 

ACCESS, DN=DUMMY34,PDN=INVELNEW. 

ACCESS, DN=DUMMY33, PDN=BL0W2DN£W. 

ACCESS, DN=DUMMY51 ,PDN=BL2DINPNEW. 
x.  COPY  DATASETS 

x.  COPY  FILE34 

REWIND, DN=DUMMY34 . 

COPYR, I = DUMMY 34 , 0=FILE34, NR . 

REWIND, DN=FILE34. 
x.  COPY  FILE33 

REWIND, DN=DUMMY33. 

COPYR , I =DUMMY33, 0=FILE33, NR . 

REWIND, DN=FILE33. 
x.  COPY  FILE51 

REWIND, DN=DUMMY51 . 

COPYR , I =DUMMY51 ,0=FILE51,NR. 

REWIND, DN=FILE51 . 
x . 

STATIC, LEVEL=NEW. 

BUILD, 1=0, L -0, OBL=KCMCAIR, B=0, NBL=MCARBIN, NODIR .  CREATE  SEQUENTIAL  BIN 

x . 

ASSIGN, DN=INPUT,A=FT05.  ALLOCATE  APPROPRIATE  FILES 
ASSIGN, DN=$0UT.A=FT06 . 

ASSIGN, DN=GEOM,A=FT29. 

ASSIGN, DN=FILE34,A=FT34. 

ASSIGN, DN=FILE33,A=FT33. 

ASSIGN, DN=FILE51,A=FT51 . 


'  data  from 
,  previous 
run 


MODE, FI=DI SABLE, 
x . 

3EGLDR,CMD=,BIN=MCARBIN« .  LINK  BINARIES 
vABD .  EXECUTE  LOAD  MODULE 

x . 

SAVE, DN=FILE34,PDN=INVELNEW.  SAVE  OUTPUT  DATASET  FOR  THE  NEXT  CAL. 
SAVE, DN=FILE33,PDN=BLOW2DNEW. 

SAVE, DN=FILE51,PDN=BL2DINPNEW. 

x  . 

EXIT. 

/EOF 

FI 5  WITH  LAMINAR  BUBBLE  WING,  MULTIPLE  ALPHAS,  M=0.6 
&A  LIFT3D=1,IBL3D=0, MAXCYC=2 , MLINE=51,NLINE=14, 

MAXANG=2, ALPHAI ( I)=7. 95, ALPHA I(2)=8. 95, AMACHC 1 )  =  0 . 60 , 
AMACH(2)=0.6,AEND 

SB  BOV2= 1 3 . 008 . I AUT0W=2 , NSYM1 =1 , L IFSEC=1 , NSTRIP=1 3, 

I XFL  AG( 1 )  =  3, NS0RCE=50 , NWAKE  =  1 , RFAREA  =  96 . 7 , REFCHD=8 . 325, SEND 


saved  for 
future 


6 . 022 

8.325 

1.0 

0.0 

610 

2.0 

10.0 

0.0 

51  .  0 

2.0 

5.0 

0.0 

6 . 022 

8.325 

1.0 

0.0 

61  .0 

2.0 

10.0 

0.0 

51 . 0 

2.0 

5.0 

0.0 

input 
data  set 


/EOF 

32.966003  13.008000  5.949856821 
32.924698  13.008000  5.950995400 
32.863266  13.008000  5.952991500 
32/82425  13,008000  5.956209200 


43.000000  0.5522600  5.600000400 
43.000000  0.1862400  5.600000400 
43.000000  0.0000000  5.600000400 


geometry 
data  set 
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SECTION  V 
TEST  CASES 


Three  test  cases  were  investigated:  (1)  F-15  with  laminar  bubble  wing, 

(2)  Advanced  Navy  Fighter  (ANF),  and  (3)  unmodified  F-15. 

5.1  Test  Case  1  -  F-15  with  Laminar  Bubble  Wing 

The  coordinates  of  the  geometry  for  test  case  1  are  given  in  model  scale 
(4.7%)  in  inches.  1084  panels  defined  this  configuration  with  600  of  them 
describing  the  wing,  as  shown  in  Pigure  26.  All  of  the  calculations  use  the 
symmetry  option  of  the  program  so  that  the  given  number  of  panels  are  for  only 
half  the  configuration. 

Inviscid  and  viscous  calculations  for  case  1  were  made  to  compare  with 
test  data  [32]  at  a  Mach  number  “0.6,  Rc  -  8.5  x  106/ft  and  a  range  of  angle 
of  attack  from  6.84°  through  15.83°.  Viscous  solutions  were  calculated  by  both 
the  two-dimensional  strip  theory  boundary  layer  and  the  quasi-three-dimensional 
boundary-layer  procedures,  and  the  location  of  transition  was  calculated  by  the 
program  using  Michel's  criterion. 


The  input  required  to  execute  the  two-dimensional  strip  theory 
boundary-layer  procedure  for  case  1  was: 


1.  A  data  set.  Unit  29,  Format  (3P10.0,  211),  containing  the  x,y,z 
coordinates  of  the  geometry,  and 


2.  A  data  set,  Unit  5,  containing  the  flags,  geometry  and  flow 
parameters,  and  boundary- layer  parameters. 


FIS  -  LAMINAR  BUBBLE  WING,  A=8.84,  M=0.6,  TRANS.  CALC. 
*A  LIFT3D=1, IBL3D=0, MAXCTC=2, MLINE=51, NLINE=14, 
MAXANG=1 , ALPHAI (1) =B . 84, AMACH=0, 60, SEND 
SB  B0V2=13 . 008, IAUT0H=2, NSVM1=1, LIFSEC=i, NSTRIP=13, 
NS0RCE=50, NWAKEsl,  RFAREA=98 .7, REFCHD=8 . 325, SEND 
6.022  8.325  0.0  0.0 

81.0  2.0  10.0 

51.0  2.0  5.0 


If  more  calculations  are  desired  for  additional,  higher  angles  of  attack, 
the  results  of  the  calculations  of  a  lower  angle  of  attack  may  be  used  for  the 
Initial  guess  for  these  calculations.  The  data  saved  on  units  33  and  51  from 
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the  first  calculation  is  now  used  as  input  for  the  additional  calculations  and 
the  namelist  block  input  will  appear  as: 


F15  WITH  LAMINAR  BUBBLE  WING,  MULTIPLE  ALPHAS,  M=0.6 
AA  LIFT3D=1, IBL3D=0, MAXCYC=2, MLINE=51, NLINE=14, 
MAXANG=2, ALPHAI (1) =10 . 75, ALPHAI (2) =12 . 95, 

AMACH=0 . 6, (END 

AB  B0V2=13 . 008, IAUT0W=2, NSYM1=1, LIFSEC=i, NSTRIP=13, 
NS0RCE=50, NWAKE=1,  RFAREA=96 . 7, REFCHD=8 . 325, AEND 


6.022 

8.325 

1.0 

Q.O 

61.0 

2.0 

10.0 

51.0 

2.0 

5.0 

6.022 

8.325 

1.0 

0.0 

61.0 

2.0 

10.0 

51.0 

2.0 

5.0 

To  execute  the  program  using  the  quasi-three-dimensional  boundary- layer 
procedure,  the  input  is  identical  to  that  of  the  two-dimensional  strip  theory 
boundary- layer  procedure  with  the  exception  of  the  following  flags  in 
Namelist  A:  IBL3D  =  1,  INTF  =  1. 

5.2  Test  Case  2  -  Advanced  Maw  Fighter  (AHF) 

Geometry  data  for  the  Advanced  Navy  Fighter  consisted  of  6%  model  scale 
coordinates  given  in  inches.  The  configuration  (Figure  27)  was  defined  by 
approximately  1200  panels  which  provided  for  360  panels  to  describe  the  wing 
and  192  panels  to  describe  the  canard.  Since  the  boundary- layer  effects  can 
be  calculated  on  only  one  lifting  section,  this  configuration  was  first  run 
without  the  canard  (Figure  27a)  and  the  results  were  compared  with  experiment 
(Figures  15  and  16).  Calculations  with  the  canard  on  (Figure  27b)  were  then 
made  and  the  effect  was  compared  with  the  canard  off  case  (Figure  17). 

Because  the  wing  is  mounted  so  low  on  the  fuselage  the  wing  extra  strip 
would  protrude  through  the  fuselage  if  it  were  defined  as  a  continuation  of 
the  wing  to  the  centerline  as  is  usually  done.  Therefore  the  extra  strip  was 
tilted  "up"  so  that  it  fell  completely  inside  the  fuselage  so  as  not  to  change 
the  fuselage  geometry.  See  Figure  28. 

Calculations  for  test  case  2  were  made  for  M  =  0.5988,  R  =  7.0xl06/ft, 

u 

C  -  0.45  and  a  =  7.7182°  to  compare  pressure  data  given  in  Reference  34. 

L 

Viscous  solutions  were  obtained  with  the  strip  theory  boundary- layer  procedure 
being  used  to  calculate  the  viscous  effects.  As  in  test  case  l,  the 


transition  location  was  calculated  by  the  program  using  Michel's  criterion. 
The  namelist  block  input  for  this  case  was: 

ANF.  NO  CANARD . , TRANS .  CALC.,  2DSTRIP  BL,  h=0.6 
(A  LIFT3D=0, IBL3D=0, 

MAXCYC=2, NLINE=8, nLINE=51, 
t1AXANG  =  l,  ALPHAI  (1)  =  7. 7 182,  At1ACH=0. 6,  SEND 
tB  B0V2=12. 955, IAUT0W=2, NSYI1i  =  l, LIFSEC=1, REFCHD=8 . 2912, 

LINEAR  =  0, IXFLAG (1) =3, IPCY  =  1 , 

RFAREA  =  95 . 904, NS0RCE=50, NSTRIP=7,  NWAKE=1,  tEND 
7.0  12.0  0.0  0.0 

70.0  2.0  10.0 

51. 0  2.0  5.0 


A  series  of  calculations  at  various  angles  of  attack  were  made  to  compare  the 
calculated  and  experimental  lift  coefficients  for  the  canard-off  case.  These 
calculations  were  repeated  for  the  configuration  with  the  canard  to  compare 
the  calculated  effect  of  the  canard.  The  namelist  block  input  for  the 
configuration  with  canard  for  multiple  angles  of  attack  was: 


ANF.  CANARD. , TRANS.  CALC.,  2DSTRIP  BL,  f1=0.G,  nULT.  ALPHAS 
tA  LIFT3D=1, IBL3D- 0, 
riAXCVC=2,NLINE=6,  MLINE  =  S1, 

rtAXANGsB,  ALPHAI  (1)  =1 . 3,  ALPHAI  (2)  =3 . 4,  ALPHAI  (3)  =4 . 7, 

ALPHAI  (4)  =6 . 0,  ALPHAI  (5)  =7 . 7182,  ALPHAI  (6)  =9 . 2,  Af1ACH=0 .6,  tEND 
tB  B0V2=12 . 955,  IAUT0W=2,  NS¥f1i=l,  LIFSEC  =  2,  REFCHD  =  8 . 2912, 
LINEAR=0, IXFLAG (1) =3, IXFLAG (2) =3, IPCV=1 , RFAREA=95 . 904, 

NS0RCE (1) =50,  NS0RCE (2) =48, NSTRIP (1) =7, NSTRIP (2) =5, NWAKE (1) =1, 
NWAKE (2) =1, tEND 


7.0 

12.0 

0.0 

0.0 

70.0 

2.0 

10.0 

51.0 

2.0 

5.0 

7.0 

12.0 

0.0 

0.0 

70.0 

2.0 

10.0 

51.0 

2.0 

'■.0 

7.0 

12. 0 

0.0 

0.0 

70.0 

2.0 

ia.  a 

51.0 

2.0 

5.0 

7.0 

12.0 

0.0 

0.0 

70.0 

2.0 

ia.  a 

51.0 

2.0 

5.0 

7.0 

12.0 

0.0 

0.0 

70.0 

2.0 

10. 0 

51.0 

2.0 

5.0 

7.0 

12.0 

0.0 

0.0 

70.0 

2.0 

10.0 

51.0 

2.0 

5.0 

5.3  Test  Case  3  -  Unmodified  P-15 


The  unmodified  F-15  geometry  coordinates  were  full  scale  in  inches.  This 
configuration  [35]  was  defined  by  981  panels,  450  of  which  described  the  wing  (F 
ig.  29).  The  wing  is  mounted  high  on  the  fuselage  and  as  in  the  case  of  the 
Advanced  Navy  fighter  configuration,  the  extra  strip  had  to  be  tilted,  in  this 
case  "down",  so  as  not  to  protrude  through  the  top  of  the  nacelle.  Figure  30 
shows  the  arrangement  of  the  panels  on  the  wing  used  in  the  calculations  of 
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test  case  3.  The  measured  data  was  obtained  with  a  sting-mounted  fuselage/wing 
configuration  which  included  the  vertical  tail.  The  computed  results,  however, 
were  for  a  configuration  without  the  vertical  tail. 


Viscous  and  inviscid  calculations  were  made  for  this  configuration  to 

6 

compare  with  experiment  [35]  at  M  =  0.6,  R  /ft  =  5.78  x  10  for  a  range  of 
angle  of  attack,  2.21°  <  a  <  8.66°.  The  two-dimensional  strip  theory 
boundary- layer  procedure  was  used  in  the  interaction  calculation  to  compute 
the  viscous  results.  The  namelist  block  input  for  test  case  3  follows: 


NEW  FIS  CONFIGURATION,  Ms0.8,  STRIP2D  BL,  TRANSITION  INPUT 
AA  LIFT3D=1,IBL3D=0, 
tlAXCYC=2, NLINE=11, hLINEsBi, 
riAXANG  =  2,  ALPHAI  (1)  =4 . 38,  ALPHAI  (2)  =6 . 51, 

Af1ACH=0. 6,  KINKN0=1,  HS  =  6,  AEND 
AB  BOV2=256 . 20, IAUT0W=2, NSVMl=i, LIFSEC=i, 

IXFLAG (1)  =3,  IPCV=1, REFCHD=191 . 23, RFAREA=43776 . , NS0RCE  =  60, 


NSTRIP 

=10,NWAKE=1, 

AEND 

4.33 

191.23 

0.0 

0.0 

71.0 

0.0 

10.0 

6.0 

254.0 

220.0 

185.0 

150.0 

115.0 

80.0 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

61.0 

0.0 

5.0 

6.0 

2S4.0 

220.0 

185.0 

150.0 

115.0 

80.0 

-0.0S 

-0.05 

-0.05 

-0.05 

-D.05 

-0.05 

4.33 

191.23 

1.0 

0.0 

71.0 

0.0 

10.0 

6.0 

254.0 

220.0 

185.0 

150.0 

115.0 

80.0 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

61.0 

0.0 

5.0 

6.0 

254.0 

220.0 

185.0 

150.0 

115.0 

80.0 

-0.05 

-0.05 

-0.05 

-0.05 

-0.05 

-0.05 

I 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


An  interactive  viscous/inviscid  procedure  has  been  developed  for  the 
computation  of  viscous  flow  over  3-D  wing/body  configurations.  The  method 
developed  makes  use  of  a  three-dimensional  surface-source  panel  method 
developed  by  Hess  [11]  and  an  inverse  finite-difference  boundary- layer 
procedure  of  [6]. 

Two  alternative  boundary- layer  formulations  have  been  considered,  the 
first  is  a  two-dimensional  strip  theory  implementation  [4]  and  the  second  [10] 
is  a  quasi- three-dimensional  boundary- layer  method.  These  procedures  have 
been  coupled  with  the  external  potential  flow  calculation  via  an  interface 
routine  which  handles  the  necessary  interpolation  and  transfer  of  data  between 
the  viscous  and  the  inviscid  parts  of  the  calculation. 

The  present  method  is  able  to  compute  flows  with  both  leading-  and 
trailing-edge  separations.  In  most  of  the  cases  considered  here,  the  calcu 
lations  were  performed  by  computing  transition  by  an  empirical  formula,  which, 
at  higher  angles  of  attack,  resulted  in  a  separated  region  as  much  as  3b% 
upstream  of  the  trailing  edge.  However,  at  the  conditions  approaching  stall, 
a  small  change  in  transition  location  led  to  regions  of  massive  flow  separa 
tion  on  the  wing,  which  the  present  method  handled  without  any  numerical 
difficulties. 

The  studies  conducted  here  indicate  that  at  lower  angles  of  attack  there 
is  very  little  difference  between  the  results  computed  by  either  of  the  two 
boundary- layer  procedures.  However,  at  higher  angles  of  attack,  for  which 
there  are  significant  regions  of  flow  separation,  the  quasi-  three-dimensional 
boundary  layer  procedure  has  been  found  to  provide  a  useful  improvement  in  the 
results  when  compared  with  the  experimental  data. 

The  inviscid  method  used  here  is  general  and  can  be  applied  to  complex 
aircraft  configurations.  The  viscous  effects,  however,  are  limited  to  the 
wing  alone  with  some  approximations.  For  example,  the  use  of  either  two- 
dimensional  strip  theory  or  the  quasi-three-dimensional  boundary- layer  approx 
imation  is  satisfactory  at  low  to  medium  angles  of  attack.  At  higher  angles 


4f. 
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of  attack,  the  flow  becomes  "more"  three  dimensional  and  the  calculations 
should  use  the  complete  equations.  Furthermore,  the  effect  of  the  viscous 
wakes  neglected  in  the  solution  procedure  also  becomes  important  at  higher 
angles  of  attack  and  needs  to  be  included  in  the  solution  procedure. 

The  viscous  capabilities  of  the  present  computer  code  can  also  be  extended 
to  the  other  components  of  the  airplane  configurations  to  include,  for  example, 
multiple  lifting  sections,  as  well  as  the  fuselage  and  other  nonlifting 
components . 
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Figure  1.  Wing/body  configuration. 


Figure  2.  Effect  of  alternate  Kutta 
conditions. 
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Figure  3.  Effect  of  wake  on  the  separation  region  and  displacement  thickness 
for  the  NACA  0012  airfoil,  Rc  ,  6.0  x  106. 
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Figure  6.  Effect  of  transition  location  on  flow  separation,  displacement 
thickness  and  lift  distribution  for  the  RAE  wing,  a  =  17.5°. 
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Figure  7.  Computed  transition  location,  flow  separation,  displacement 
thickness  and  lift  distribution  for  the  RAE  wing,  u  =  18.5° 
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Figure  11.  Effect  of  transition  location  on  flow  separation,  displacement 
thickness  and  lift  distribution  for  the  NACA  0012  swept  wing, 
a  =  21.12°. 


Figure  12.  Comparison  of  2-0  strip  theory  and  quasi-3-D  methods  for  the  F-15 
fighter  with  laminar  bubble  wing. 
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Figure  13.  Comparison  of  the  calculated  and  experimental  pressure  distribution 
for  the  F-15  with  laminar  bubble  wing,  a  =  10.75°.  Insert  shows 
the  chordwise  variation  of  displacement  thickness  as  calculated  by 
the  2-D  strip  theory  and  the  quasi -three-dimensional  boundary- 
layer  methods,  (a)  38%  semi  span,  (b)  58%  semi  span. 
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Figure  13.  (Continued)  (c)  88%  semispan. 
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Figure  14.  Comparison  of  the  mid-semispan  pressure  distribution  obtained  by 
the  panel  method  and  the  transonic  flow  method  for  the  F-15  with 
laminar  bubble  wing,  M  =  0.6.  (a)  a  =  6.84°. 


Figure  15.  Comparison  of  computed  and  experimental  lift  curves  for  the  ANF 
wing/body  configuration. 
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Figure  19. 


Mgure  ?0.  Flow  structure  of  three-dimensional  viscous/inviscid  interaction 
program. 
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Figure  21.  Examples  of  correct  and  incorrect  input. 
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Plan  view  of  the  input  points  on  a  body  divided  into  sections. 
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Figure  25.  (Continued)  Output  from  the  sample  case. 
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Figure  25.  (Continued)  Output  from  the  sample  case. 
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Figure  25.  (Continued)  OutPu 
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Figure  25.  (Continued)  Output  from  the  sanple 
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Figure  25.  (Continued)  Output  from  the  sample  case. 
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Figure  25.  (Continued)  Output  from  the  sample  case. 
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Figure  25.  (Continued)  Output  from  the  sample  case. 
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Figure  25.  (Concluded)  Output  from  the  sample  case. 


